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NONLINEAR ANALYSIS OF BONDED JOINTS 
WITH THERMAL EFFECTS 


ABSTRACT 

A numerical analysis of the nonlinear response of bonded joints is 
presented. Mechanical and thermal loadings are considered. Material 
stress-strain response is represented by Ramberg-Osgood approximations. 
Temperature dependent properties including modulus percent retentions 
and coefficients of expansion are modeled with linearly segmented 
curves. Bonded Joints with graphite-polyimide, boron-epoxy, titanium, 
or aluminum adherends are analyzed using a quasi 3-dimensional finite 

element analysis. In adhesive bonded Joints, the adhesives considered 
are Metlbond 1113 and AF-126-2. 

Elastic results are presented for single and double lap joints, with 
and without adhesives. It Is shown that mechanically induced stresses 
are greatly affected by longitudinal adherend stiffness. The effects 
of adherend transverse stiffness are shown to be significant in some 
cases. Residual curing stresses are shown to be significant in all 
joints except those with similar adherends and no adhesive. 

Nonlinear results are presented for adhesive bonded Joints. It is 
shown that adhesive nonlinearities are only significant in the predicted 
adhesive shear stresses. Adherend nonlinearities and temperature de- 
pendent properties are shown to have little effect upon the adhesive 
stress predictions under mechanical and thermal loadings. 
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Chapter 1 
INTRODUCTION 

With the development of advanced fiber reinforced composites as 
viable structural materials the adhesive bonded joint has become of 
primary importance. The bonded joint does not require that the stru- 
ctural members being joined (adherends) be perforated to facilitate 
bolts. Without the bolt holes and the stress concentrations associated 
with them, a substantial weight savings can be realized which is a major 
reason for selecting composite materials for the structural component. 

In order to fully realize the strength of the composite adherends, 
the adhesive bonded joint must be efficiently designed and this requires 
an adequate prediction of the stress distributions in the adhesive layer 
of the joint. The study of stresses in the adhesive layer has been 
approached by researchers in the past using one of two types of analysis 
procedures. 

Many researchers have attempted to predict the stresses using a 
closed form analytical solution. However, when using this approach, the 
equations that need to be solved become exceedingly complicated and this 
leads to the need for simplifying assumptions. These assumptions have 
included linearity, material isotropy, restrictions on the geometry of 
the joint, and neglect of thermal effects. 

Other researchers have approached the problem through the use of 
numerical techniques such as finite element analysis. They have usually 
found it necessary to use some or all of the assumptions made for the 
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closed form type solution procedure. The motivation for the present 

study is to show the capability of a finite element computer analysis 

program developed by Renieri and Herakovich [1] to adequately predict 

these stress fields. The program has the capability for material ortho 

tropy, material nonlinearities, and temperature dependent properties. 

Modifications to the computer program for this study include increased 

element capacity, improved execution time, and capability for hygrother 
mal analysis. 
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Chapter 2 
LITERATURE SURVEY 

The first investigation into the behavior of bonded joints found in 
the literature was presented in 1944 by Goland and Reissner [2]. They 
obtained an analytical solution by assuming a state of plane strain, 
prescribing the distributions of the shear and peel stresses to be 
parabolic and linear, respectively, and applying restrictions on the 
ratios of adherend moduli and thickness to adhesive modulus and thick- 
ness. The solution is based on the principle of minimum potential 
energy and is restricted to linear isotropic materials and identical 
adherends. Thermal effects were neglected for this analysis. 

Erdogan and Ratwani [3] approached the problem of an orthotropic 
plate bonded to an isotropic plate with an isotropic adhesive. They 
obtained closed fom solutions for stepped lap and smoothly tapered 
joints. Their solution was based on a summation of forces in the ad- 
hesive layer and adherends and assumed plane stress and linear material 
behavior. This solution predicted stress singularities at the edges of 
each step in the adhesive in the stepped lap joint, and at the ends of 
the overlap in the adhesive for the smoothly tapered joint. 

Barker and Hatt [4] used a linear elastic finite element computer 
analysis program to compare results with the work of Erdogan and Ratwani 
[3]. The adherends were modeled using four noded isoparametric elements 
and material orthotropy was considered. The adhesive layer was modeled 
using a special element developed for that purpose which was formulated 
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to have no thickness and used modified stiffnesses derived from the 

moduli, thickness, and length of the adhesive layer. Their results ^ | 

compared favorably with Erdogan and Ratwani. I 

Sainsbury-Carter [5] solved the stepped and linearly tapered “ | 

bonded joints by assuming linear isotropic materials and solving the \ 

equations of equilibrium. It was assumed that the moduli of the adher- i 

ends are much larger than those of the adhesive. Also the analysis was ' 

one-dimensional and thermal effects were neglected. It was shown that i 

j 

the thickness of the adherends greatly affected the magnitude of the I 

> 

peak shear and peel stresses and an iterative technique v-zas developed to I 

modify these thicknesses within stress design criteria. 

i 

Wah [6] investigated non-symmetric single lap joints with composite ; 

adherends and isotropic adhesives. Laminated plate theory was used to 
develop stress and moment resultants and relate them to mid-plane ,i 

strains and curvatures. The laminated adherends were required to be I 

mid-plane syirnietric in order to uncouple the bending-stretching terms in ^ 

the previously mentioned relationships. A solution was presented for i 

the joints under shear loadings as well as axial loads. The solutions ; 

were restricted to the elastic range and thermal effects were neglected. ^ 

Hart-Smith [7] was the first researcher to consider the non-lin- 
earity of the adhesive layer by assuming its stress-strain response to ; 

be elastic-perfectly plastic. This effort presents solutions and ^ I 

design aids, including thermal effects, for single, double and stepped 

i: 

lap joints and linear tapered joints. While a thermal mismatch is "" i 

j 

considered all material properties are considered temperature indepen- I 
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dent. This research notes the increased failure strength by allowing 
for plastic deformation in the adhesive. 

A comparison of theoretical and experimental shear stress in a 
single lap joint was presented by Sharpe and Muha [8]. The joint 
modeled had plexiglass adherends and an epoxy adhesive. Good correla- 
tion was obtained with the work of Goland and Reissner [2] and with a 
linear computer analysis program, B0ND4, of the University of Delaware. 

Renton and Vinson [9-10] performed parametric studies on single lap 
joints as well as fatigue testing and thick adherend lap shear testing. 
The specimens were comprised of mid-plane symmetric composite adherends 
and elastic behavior only was studied. Linear thermal effects were also 
included. The parameters studied were over-lap length, adhesive thick- 
ness, and ply orientation in the composite adherends. Comparisons were 
made with the work of Goland and Reissner [2] with Renton and Vinson's 
work showing better satisfaction of stress free boundary condition at 
the edge of the adhesive layer. 

Grimes, Greimann et al [11] approached the analysis of single, 
double, and stepped lap joints from both the finite element method and 
numerical integration of the governing differential equations. Their 
analysis included full material nonlinearity in the adherends and 
adhesive layer. The development for both solutions was based on the 
deformation theory of plasticity with the finite element analysis 
utilizing an iterative procedure until the solution converged. In both 
forms of analysis, solutions were presented for room temperature only 
and curing stresses were neglected. 
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DasGupta and Sharma [12] used an analysis similar to Goland and 
Reissner [2] to predict stresses 1n lap joints with prebent adherends. 
The work showed a decrease In peak stresses with the use of bent ad- 
herends. 

Renton [13] provided an analysis of the thick adherend lap shear 
test using the work of Renton and Vinson [9-10]. This research verified 
the validity of the test. 

Other researchers have investigated the effects of moisture [14], 
and the reliability [15] of lap joints. 

While this survey is by no means all-inclusive, it is representa- 
tive of the research that has been performed and from this survey the 
need for fully-noni inear material behavior and temperature dependent 
properties can be seen as these physical realities have been consis- 
tently neglected. 


Chapter 3 

BASIC CONSIDERATIONS 

The bonded joints selected for analysis 1n this study are single 
and double lap joints with and without adhesive layers. Typical geo- 
metries for joints with adhesive layers are shown in Figure 1. Composite 
and isotropic adherends are considered with nonlinear material properties 
and thermal stresses as well as temperature dependent properties. 
Hygroscopic analysis capabilities are presented but no results are 
included in this investigation due to a lack of complete consistent 
data. 

3.1 Geometric Restrictions 

For the present study it is assumed that the joint is in a state of 
plane strain (i.e. e = 0» or e = const). This is a valid assumption 
if the X dimension of the joint (Fig. 1) is large and the cross-section 
under consideration is some distance removed from contraints that are 
dependent upon the x coordinate. 

The analysis is also restricted to balanced, mid-plane symmetric 
composite adherends and laminate material properties are used for these 
components. 

3.2 Quasi 3-Dimensional Analysis 

The analysis of reference [1] considers a long prismatic bar under 
the influence of a uniform applied strain or temperature change to have 
strains independent of the x coordinate. With this assumption, the 
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strain displacement relations can be written as 
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where u, v, and w are x, y, and z displacements, reioectively, and 
through fg are unknown functions of y and z coordinates only. With the 
use of suitable mathematical manipulation, (3.1 ) can be integrated yield- 
ing the following form for the displacement fields 


u(x,y,z) = x(C^y+C 2 Z+C 3 ) + U(y,z) 

2 

v(x,y,z) = x{C^z+Cg) - |-+ V(y,z) (3.2) 

2 

w(x,y,z) = x(-C^y+Cg) - ^ + W(y,z) 

where through Cg are unknown constants and U, V, and W are unknown 
functions of y and z only. With this assumption, and neglecting body 
forces, the equilibrium equations can be written as 


^-^xz 
9y 9z 


= 0 


3a 9 t 

0 

8y 3z 




9y 


az 


= 0 


(3.3) 


0 


/ 


3.3 Bonded Joints 

The analysis developed in [1] can be used as two different formula- 
tions for the solution of bonded joints. The first of these corresponds 
to a classical plane strain solution in which = 0. This 

type of solution is entirely 2-dimensional and is the type of plane 
strain analysis used in a majority of joint analyses previously in the 
1 iterature. 

The second formulation is a more general plane strain solution 

where z f 0 but is equal to some constant For this procedure the 

components of strain y and y are assumed to be zero. The principal 

xz xy 

difference between the two formulations is that the second accounts for 
the transverse stiffness (E^) of the adherends and adhesive while the 
first does not. This second formulation is used for this analysis and 
comparisons between results obtained from the two solutions are pre- 
sented in Chapter 5. In the following sections and chapters the first 
formulation is referred to as a 2-dimensional formulation while the 
second is referred to as a quasi 3-dimensional analysis as it corresponds 
closely to the analysis of section 3.2. 

3.3.1 2-Dimensional Joint Formulation 

In the classic plane strain solution the displacements, strains, 
and therefore stresses, are independent of the x-coordinate (Fig. 1). 
Under these assumptions the displacement fields (Equ. 3.2) reduce to 


u = 0 
V = V(y,z) 
w = W(y,z) 


(3.4) 
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The nonzero strain components (e^, and have the same definitions 
as in Equ's. (3.1). 

3.3.2 Quasi 3-Dimensional Joint Formulation 

If, in a bonded joint, it can be assumed that is a nonzero con- 

stant and that Yv-, ®nd y are zero, then the displacement fields 
az xy 

(Equ's. 3.2) reduce to 

U = 

V = V(y,z) (3.5) 

w = W(y,z) 

In these equations the only nonzero constant from Equ's. (3.2), C^, has 
been renamed £ and corresponds to the uniform normal strain e„. The 
remaining strain components (g^, c^, and y^^) again have the same form 
as in Equ's. (3.1 ) . 

The assumption that the strain component normal to the plane of the 

analysis is constant is strictly valid for the case of single lap 

joints vn*th identical orthotropic adherends. It is also valid for 

double lap joints where the outer adherends are identical and exhibit 

orthotropy. The assumptions that y^, y = 0 are valid if the 

xz 

adherends are orthotropic. 

The restriction that all adherends be orthotropic is satisfied 
by all joints analyzed in this study. It should be noted however 
that two of the joints analyzed do not satisfy the first condi- 
tion as they are single lap joints with differing adherends. 
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For these joints, the solution neglects the effects of bending out of 
the plane of the joint and is similar to a membrane solution in this 
respect. 

3,4 Materia! Properties 

In order to adequately model a layered adherend layer by layer an 
excessive number of finite elements would be needed. Because of this, 
it is necessary to use laminate material properties and consider the 
composite adherends to be homogeneous orthotropic materials for the 
joint analysis. 

Obtaining laminate material properties from the literature proved 
to be impossible thus making it necessary to generate these properties 
analytically. For this generation of properties two different approaches 
were used. The stress-strain response of a laminate was predicted 
following the work of Renieri and Herakovich [1], while thermal proper- 
ties were predicted using classical lamination theory, 

3,4.1 Prediction of Laminate Stress-Strain Response 

The details of the analysis of ref. [1] will not be presented as to 
do so would be overly repetitious; however a brief outline will be pre- 
sented for completeness. 

The analysis utilizes the displacement fields of Equ's. (3.2). 
Because the laminates in question are balanced and midplane symmetric 
the analysis can be reduced to the quarter section shown in Figure (2b) 
with certain symmetry and anti-symmetry conditions. The displacement 


l 




r 



fields (Equ’s. 3,2) reduce to 


u = + U(y,z) 

V = V(y,z) (3.6) 

w = W(y,z) 

and again the constant has been renamed and corresponds to a uni- 
form applied strain. The displacement fields (3.6) along with the stress- 
free boundary conditions along the free edges, top and bottom surfaces 
and certain restrictions imposed upon the displacements by the symmetry 
and anti-symmetry conditions mentioned previously represent the boundary 
value problem to be solved by the finite element analysis. With this 
analysis and the nonlinear finite element program, the moduli and 
E can be predicted as functions of strain level. 

J" V 

3.4.2 Prediction of Laminate Thermal Properties 

Laminate thermal properties including coefficients of thermal ex- 
pansion and moduli as functions of temperature are predicted using 
lamination theory and unidirectional material properties as functions 
of temperature. Lamination theory as presented here cannot directly 
predict temperature dependent laminate properties, however, if the uni- 
directional properties used as input correspond to an elevated tempera- 
ture, then the laminate properties generated will also correspond to 
this temperature. Therefore, laminate properties can be predicted at 
discrete temperatures corresponding to the input data. 

The constitutive relations for a single, orthotropic lamina in the 
principal material coordinates are 


where the principal coordinates are shown in Fig. 2a as the 1-2-3 system. 
These equations can be written in an abbreviated form as 

io}^ = [C]{e}^ (3.8) 

For a coordinate rotation about the 3 axis through an angle e (Fig. 2a) 
the stresses and strains are transformed according to the following 
relations, 


= [T^]{a>^ and (3,9) 

where 



0 0 0 (m -n ) 

and m = cose, n = sine 

Combining Equ's. (3.8) and (3.9) yields 

{a}^ = [T^][C][T2]-^{e}j^ 

or 

fa}^ = [C]{e), 

where [C] is defined as 

[C] = [T^lCClCTg]"^ 


and has the form 
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‘ -1 


Equ's. (3.10) are the constitutive relations for an orthotropic lamina 
rotated through an angle e. 

For lamination theory it is assumed that a lamina is in a state of 
plane stress. It should be noted that this lamination theory development 
is the only case in which plane stress will apply while plane strain is 
assumed in all other developments. 

The mathematical statement of the plane stress assumption is 


a 


z 



e 0 


(3.13) 


which can be used to reduce the constitutive relations (Equ's. 3.7) to 



(3.14) 


This simplified form of the stiffness matrix [C] is known as the reduced 
stiffness matrix [Q]. 


If the transformation matrices (Equ's. 3,9) are reduced similarly* a 
rotated plane stress constitutive relation [Q] can be formed in the 
same manner as Equ's. (3.10). Thus 


where 


(3.15) 





^12 

^16 

■^12 

Q22 

^26 

^16 

^26 

^66 




Now, taking the standard plate theory assumption that normals to 
the mid-plane of the plate remain normal after loading, the strains of 
Equ's. (3.15) may be written as 


{e}^ = it^} + Z{ic} 

A 


where 


{e}^ = total strains 

X 


it°} = /e° \ = mid-plane strains 

*7 


y 


xy. 


(k) = / K \ = plate curvatures 


‘xy. 


and z = distance from the mid-plane 
Defining stress resultants 


{N> = / N 



xy ; \ xy 

and combining Equ's. (3.15), (3.16), and (3.17) yields 

H _ TH „ 


{N> 


= f ^ [Q]{e‘>}dz + f [q]z{K}dz 


( 3.1 


( 3.1 


( 3.1 



{N} - [A]{e°} + [B](k} 

( 

where 


[A] = 1 J 

k=1 ^ 

and 

( 

™-rJ, 


where the k denotes the number of the ply and the h's 
2b. 

are as in Fig 

For symmetric lamiantes [B] = 0 and 


{N> = [A]{e°} 

(: 

Inverting this relationship yields 


{£“} = [A]"''{N} 

C 

Noting that 


{a}^ = {N} 

{; 

and combining Equ's. (3.22) and (3.23), leads to 


{e“} = [a*]{5l ■ 

{: 



[a*J = 2H[A]“' 


These relations can be used to define the laminate properties 


E 

" ^ ^1 


E =^ = i- 

y a* 
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(3.2 


xy 




'xy 


% 


V = 4=^2 

xy 


For an orthotropic material the coefficients of thermal expanslor 


are 



(3.2i 


These coefficients transform, under the rotation defined earlier, In 
the same manner as the strains (Equ. 3,9). 



For the 2-dimensional analysis Equ. (3.27) reduces to 


a 


X 


a. 


a 


xy 


, 2 
n + m 

2mn(a^-a2) 


(3.2fi 


Now define laminate coefficients of thermal expansion such that 


{e°} = {a}AT (3.29) 

where AT is a uniform temperature change and 



Combining Equ's. (3. IS), (3.22), and (3.29), and considering symmetric 
laminates only yields 




{SDaT = [A]‘^ 



H . 

[Q]{ct}dzaT 


(3.30 


or 

{5} = [A]’’’ s CQ]''{a}'^{h^-h|^_.,) (3.31 

k“l 

When the moduli and thermal coefficients used to calculate the [C] 
matrix (Equ's. 3,7) are those corresponding to an elevated temperature, 
then the moduli of Equ's. (3.25) and coefficients of expansion of Equ's 
(3.31) will be laminate properties also corresponding to that tempera- 



f 
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ture. These reduced moduli will be used as Input in the finite element 
program. 

These methods of generating material properties were resorted to 
because of the lack of consistent experimental data found in the 
literature. 

3.5 Finite Element Formulation 

As in section 3.2.1 the presentation of the complete formulation 
would be a duplication of the work of Renieri and Herakovich [1], 
therefore only the highlights will be given here. 

The finite element solution process involves the subdivision of a 
structure into a finite number of smaller elements (Fig. 3). This 
process is known as discretization. For each of these finite elements 
a set of interpolation functions are chosen to represent the displace- 
ments at any point in the element as functions of the displacements at 
the corners or nodes of the element. Using the strain-displacement 
relations (Equ. 3.1) the strains can also be calculated as functions of 
nodal displacements. Now with the use of a variational principle, such 
as the principle of minimum potential energy, a set of equations re- 
lating nodal forces to nodal displacements can be obtained for each 
element, 

(3.32) 

where {F> = nodal forces 

(ul = nodal displacements 
[K] = element stiffness matrix 













and the superscript (£) refers to the individual element. These elemental 
relations (Equ's. 3.32) are then combined or assembled into a larger 
system of equations relating forces to displacements for the entire 
structure. The solution of these equations yields the displacements and 
therefore the strains and stresses over the entire body. 

The finite element scheme developed by Renieri and Herakovich [1] 
utilizes the constant-stress, constant-strain triangular element with 
three nodes. The interpolation functions used are 

u = a^ + a 2 y + agZ + 

v = + agy + agZ (3.33) 

w = ay + agy + agZ 

As can be seen from the form of Equ's. (3.33) these are linear relations 
and will yield constant strains when substituted into Equ's. (3.1). The 
constants a-j through 3g are functions of the spatial coordinates and 
nodal displacements of the individual elements and ? is the applied 
uniform strain. 

Manipulation of Equ's. (3.33) and substitution into Equ. (3.1) 
yields the following strain-displacement relations for an element. 


av^ + CV2 + evg 

bw-j + dWg + gWg 

bv-j + dvg + gvg + aw^ + cwg + ew^ 

bu^ + dUg + gug 

au-j + cug + eug 

where A = area of the element 

U-] 9 Ugs = x-displacements at nodes 1, 2, and 3 respectively 

Via Vgs Vg = y-displacements at nodes 1, 2, and 3 respectively 

w-| a Wg, Wg =■ 2 -displacements at nodes 1, 2, and 3 respectively 

and a, b, c, dj e, g are known constants involving spatial coordinates 

only. 

For the case of a uniform thermal load the strains are 

{e}., = {e*}^ - (3.35) 

where {e*}^ = total strain 

{e}^ = mechanical strain 
and {e^}^ = thermal strain 

which consists of {a}^ (Equ. 3.26) multiplied by the temperature change 
AT. Transforming Equ's. (3.35) to an arbitrary coordinate system yields 
an individual element. 

= {e*}^ - [TglcJ}-, (3.36) 














Noting Equ's. (3.34), these strains can be written as 


- (m^a^ + 

(av^ + cvg + ev3)/A - (n^a^ + m^ag)AT 
(bw, + dw« + gw«)/A - h-aT 

, ' V (3.37) 

(bv-j + dvg + gv^ + aw^ + cWg + 

(bu^ + dUg + gu2)/A 

(au^ + cug + eu^)/J\ + 2mnAT(a^-a2) 

The preceding formulation for thermal strains is completely analogous 
for that of hygroscopic strains. For an orthotropic material the coef- 
ficients of hygrosopic expansion are 

/ 



(6)^ = 



(3.38) 


Following exactly the development of Equ's. (3.35) through (3.37) and 
substituting for {a}^ and aM for AT where AM is a uniform percent 
weight change due to moisture absorption or desorbtion, yields 
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5 - (m 3 , + n 3 «)aM 
X ^ 2 o 2 

(av^ + CV2 + - (n 3^ + ni 3 -|)AM 

(bw-j + dWg + gw2)/A - 33AM 

(bv^ + dV2 + gV3 + aw^ + CW 2 + eW3)/A 

(bu^ + dU2 + gu3)/A 

(au^ + CU 2 + eu3)/A + 2mnAM(3^-32^ 

The hygroscopic expressions are presented here because the capability 
for this type of analysis has been included In the computer program. 
Results will not be presented because of the lack of data as stated 
earlier. It should also be noted that the derivation is for a uniform 
temperature or moisture change and that analysis should be limited to 

cases where uniformity is a valid assumption. 

The principle of minimum potential energy states that a body is in 
equilibrium when the total potential energy is minimum where 

= Ue + Wg, (3*^0) 




U = internal strain energy 
e 


and 


Wg = potential energy of the applied loads 
The internal strain energy for an element is 

Ug = iJie)''’[C3{e}dVol 


(3.41) 


which for an element with constant strains and unit thickness reduces 



28 


U = §- {e}'''[c]{e} (3-42) 

In both Equ's. (3.41) and (3.42) the strains Involved depend upon whether 
the loading is mechanical (Equ's. 3.34), thermal (Equ's. 3.37) or 
hygroscopic (Equ's. 3.39). The potential energy of the applied loads is 
given by the negative of the applied forces multiplied by their respective 

displacements. 

Minimization of Equ's. (3.40) with respect to displacements yield 
the elemental stiffness matrix plus strain, thermal and hygrosopic re- 
lated vectors. The forms for these can be found in Appendix A. 

3.6 Boundary Conditions for Joints 

The boundary conditions applied to single and double lap joints for 
the present study are shown in Fig. 3. A number of different boundary 
conditions and loadings were investigated and comparisons were made. 

These are summarized in Table 1- Noting Table 1 it is seen that, 
all conditions except the fourth yield comparable results. The lower 
peak stresses for this condition can be attributed to an overly flexible 
model. It was reasoned that this dees not correspond to physical reality 
as a real joint would not be free to deflect up and down where the load- 
ing is applied. The first set of conditions was eliminated from consider 
ation as they can only be applied to symmetric single lap joints. The 
second set was eliminated because during the solution process a negative 
diagonal in the global stiffness matrix was encountered. The third set 
of conditions was disregarded because the stress distribution in the 
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TABLE 1 


Comparison of Various Boundary Conditions for Symmetric Single Lap Joints 








Stable 

Yes 


Symmetric Stress 
Distribution 

Applied 

Load 

Max 

■^yz 

Max 

^y 

Max 

Max End 
Displacement 

Yes 

300 

(lbs) 

2.39 

(KSI) 

.84 

(KSI) 

1.13 : 
(KSI) 

.91 X 10*3 
.82 X 10-3 
(in) 

Yes 

300 

(lbs) 

2.39 

(KSI) 

.84 

! 

1.13 

.91 X 10-3 
.82 X 10-3 

No 

1 

! 

1 

300 

(lbs) 

j 

'2.52 

2.25 

(KSI) 

.97 

.72 

(KSI) 

1.32 

.93 

(KSI) 

.179 X 10“2 
.165 X 10-2 
(in) 

Yes 

.165 X 10"^ 
(in) 

1.91 

.66 

.88 

.165 X 10"^ 
.165 X 10-2 

Yes 

i 

.165 X 10-2 
(in) 

2.32 

(KSI) 

.83 

(KSI) 

1.11 

(KSI) 

.165 X 10“^ 
.165 X 10-2 
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adhesive layer was not symmetric for a symmetric joint while it is in- 
tuitively obvious that it should be. The fifth set of constraints and 
loading was chosen because it provided symmetric results, numerical 
stability and is applicable to nonsyrametric joints. 

3.7 Qualifying Notes 

Due to the lack of data present in the literature a number of 
assumptions concerning material properties have been made. For lamina 
data input, it is assumed that Egg is identical to E33. It is also 
assumed that the three shear moduli G^g, 6^3 and Ggg are identical. 
Poisson ratios are considered constant and percent modulus retentions 
are assumed identical in tension and compression 

When inputting laminate material properties, much of the data is 
generated according to the analyses presented earlier. It is assumed 
that the Isniinate response is identical to the Og-cg response of 

a unidirectional lamina. It is also assumed that the 3^^ 

^xz"^xz ^^^ves are the same as the ti 2“"^12 of a lamina and 

that the laminate poisson ratios and v are the same as the lamina 
poissons ratios v-jg and Vgg, respectively. Other restrictions are the 
same as for lamina data. 

It should be mentioned here that the material properties used in 
this analysis {Appendix C) are not consistent. The data has been taken 
from a number of sources and is not all realted to identical material 
systems. Even with these limitations it is felt that the investigation 
still fulfills its goal of showing the capability to analyze joints if 
consistent properties were available. 








Chapter 4 

NONLINEAR ANALYSIS 

A significant portion of the research effort involved in this study 
was devoted to the modification of the finite element analysis program 
NONCOM [1]. The improved version, NONCOMl will be briefly described 
here. 


4.1 Modifications of NONCOMl 

For the analysis of bonded joints, a number of modifications to the 
existing finite element program were implemented. These included: 

(1) Increased finite element capacity 

(2) Inclusion of a more efficient equation solver 

(3) Capability for input of fully 3-dimensional orthotropic 
material properties 

Other modifications included: 

(4) Capability for hygroscopic analysis 

(5) Capability for elevated temperature and moisture 
content analysis 

4.1.1 Increased Finite Element Capacity 

In order to model an adhesive bonded joint a large number of finite 
elements are needed because of the large stress gradients and inherent 
large aspect ratio of the adhesive layer. For this reason the maximum 
number of elements was increased from 100 to 400 elements. This was 
done with an increase of high speed storage of approximately 50 percent. 
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This increase was held this low through the use of storage addressing 
schemes where many arrays are stored 1n the same common storage loca- 
tions and by the block storage schemes of the equation solver chosen. 

4.1.2 Solution of Simultaneous Equations 

The equation solver in NONCOM is .lot suitable for large systems of 
equations as its solution time becomes excessively large for such 
systems. For this reason a new equation solver, SESOL [16] was selected 
for NONCOMl . This equation solver offers a fast solution time and high 
speed storage reduction schemes. In brief the solution algorithm 
considers the system of linear equations 

[K]{X) = {R} (4.1) 

where [K] is the assembled stiffness matrix, {X} is the nodal displace- 
ment vector and {R} is the applied nodal load vector. The stiffness 
matrix is factored into an upper and lower triangular matrix 

[K] = [if [G] (4.2) 

where [G] is upper triangular and [L]"^ is lower triangular and normalized 
such that L.^. = 1 (i not summed)- Since [K] is symmetric 

G— = summed) (4.3) 

and equation (4.2) can be written as 


[K] = [L]T[D][L] 


(4.4) 
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where 

[D] = G.. (i not summed) (4,5) 

Defining 

{V} = [D][L]{x} (4.6) 

and combining equations (4.1) and (4.6) yields 

[L]'^{v} - {R} (4.7) 

The vector {v} in equation (4.7) is first found by Gauss reduction of 
the load vector and then the nodal displacements are found through back 
substitution into equation (4.6), The stiffness matrix and load vector 
are assembled and stored on low speed storage in block form as in Fig. 
4. During the solution process, reductions are performed on non-zero 
terms only and only two blocks need be in high speed storgage at any 
time. 

In the computer program NONCOMl the number of equations per block 
is determined as a function of the maximum half bandwidth plus the 
diagonal of the assembled stiffness matrix. This is done to minimize 
the number of blocks necessary and maximize the number of equations per 
block within high-speed storage limitations. By making the number of 
blocks a minimum, I-O operations performed by the computer with the 
elemental stiffness matrices are also minimized. To further this re- 
duction of operations the blocks are assembled two at a time. 
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Block Storage of Stiffness Matrix and Load Vector 
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It should be noted the size of the maximum half bandwidth plus 
diagonal is not only a function of node numbering, as is usually the 
case with the finite element analysis, but is also a function of the 
type of loading applied. When the applied loading is thermal, hygro- 
scopic or an average force applied in the x direction an equation re- 
lating a uniform strain to an average force in the x direction is 
required in addition to the equations relating nodal forces and dis- 
placements. This equation is related to all of the elements and may 
involve all the nodal displacements of the finite element model. To 
minimize the effect of this equation on the maximum half bandwidth plus 
diagonal it is assembled in the center of the stiffness matrix. The 
average force equation and its effect on the bandwidth can be seen in 
Fig. 4 where the T's represent the average force terms. The half band- 
width plus diagonal without the average force equation is shown by the 
dotted line. These two bandwidths are for identical finite element 
models under different loadings. 

4.1.3 Three Dimensional Properties 

The program NONCOMl was given the capability for fully three 
dimensional orthotropic material properties because of the restriction 
that composite adherends must be modeled as homogeneous orthotropic 
laminates (Chapter 3). This means that a plane of transverse isotropy 
cannot in general be assumed as done in NONCOM [1]. 

4.1.4 Hygroscopic Analysis 

A capability for hygroscopic loadings has been included in NONCOMl. 
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It is modeled in the computer program in exactly the same fashion as the 
thermal analysis and is therefore subject to the same limitations. The 
most notable of these is the restriction that the moisture distribution 
be uniform throughout the finite element model. Other restrictions will 
be mentioned in the section dealing with the nonlinear analysis. 

4.1.5 Elevated Temperature and Moisture Content Analysis 

This modification allows for mechanical loading at elevated tempera- 
ture and moisture content or thermal loading at elevated moisture 
content or hygroscopic loading at elevated temperature. Insight into 
the interactions between thermal and hygroscopic material response could 
not be obtained from the literature. For this reason two assumptions 
are made about these interactions. First, it is assumed that hygro- 
scopic properties are independent of temperature and that thermal proper- 
ties are independent of moisture content. It is also assumed that 
changes in mechanical properties due to temperature and moisture content 
are cumulative. By cumulative it is meant that when a modulus is to be 
modified to correspond to both temperature and moisture content it is 
first modified for the temperature and then this new modulus is then 
modified to correspond to the moisture content. The process of changing 
material properties will be more fully described in the following section. 

4.2 Nonlinear Analysis 

In order to simulate nonlinear material behavior in a computer 
program with the linear elastic finite element analysis described in 
Chapter 3 two separate' problems must be dealt with. First a method of 
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accounting for material nonlinearities within the analysis procedure and 
second a method of representing the nonlinear material properties. 

4.2.1 Incremental Loading Procedure 

The finite element computer program NONCOMl deals with varying 
material properties through the use of a Incremental solution procedure. 
With this type of procedure, the load whether mechanical, thermal, or 
hygroscopic is applied as a series of increments. This yields a series 
of linear solutions with total stresses, strains, and displacements 
formed by summation of the linear increments of these quantities. When 
applying the load incrementally the material properties are updated to 
correspond to the current levels of strain, temperature, and mositure 
content. With the finite element method an individual element can have 
material properties varying independently of other elements in the 
model . 


4.2.2 Nonlinear Data Input 

Material stress-strain response for an othrotropic material in the 
principal material coordinates are represented in the form of modified 
Ramberg-Osgood [17] approximations which have the form 

E == f + K-a”i 1 = 1 or 2 (4.8) 

In equation (4.8) E is the elastic modulus and K. and n^. are Ramberg- 
Osgood coefficients. A method for calculating the four coefficients K^. 
and n. is described in Ref, [1]. A tangent modulus can be defined as 
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p_ da 
de 


K.En.a"rVl 


i = 1 or 2 


(4.9) 


where I is the tangent modulus corresponding to the principal stress a 

Noting F,g. (5) the value of the stress corresponding to the strain 
at the end of load increment P is 


P . 
k=l 




(4.10) 


where Ae^ is the increment of strain during the load increment. 
Combining equations (4.9) and (4.10) yields for the P + l^h increment 


tP+1 - E 

p — r: — r~ i i or 2 ( 4 . 11 ) 

K.En.[ E AeJE^]"-Vi 
J=1 

With equation (4.11) and principal material strains, the tangent moduli 
are calculated at the end of each increment to be used for the next 
increment. Moduli determined are E 33 . 533 . 5 ^ 3 , and 

It can be seen in Fig. 5 that the strain for which the tangent 
modulus is calculated differs from the strain where the modulus 
should be calculated. This difference is a function of the size of the 
load increment and can be made negligible by choosing an appropriately 
small increment. For the computer analysis it is assumed that the shear 
response is independent of sign while extensional behavior can be 
different in tension and compression. 

Temperature and moisture dependent properties are represented as 
linearly segmented curves. These properties consist of percent modulus 
retention curves which represent the change in stiffness of a material 


due to variations in temperature or moisture levels and coefficients of 
thermal and hygroscopic expansion as functions of temperature and 
moisture content respectively. During thermal or hygroscopic loading, 
the moduli and coefficients of expansion are calculated at the mid-point 
of the increment using linear interpolation. For input to the computer 
program it is assumed that percent modulus retentions are identical for 
tension and compression. 


Chapter 5 

RESULTS AND DISCUSSIONS 


The analysis procedures presented in Chapters 3 and 4 were used to 
generate laminate properties and analyze various bonded joints. The 
joints investigated included single and double lap joints with, and 
without adhesive layers. Both elastic and non-linear results are pre- 
sented. The material systems considered were graphite-polyimide, boron- 

epoxy, titanium, and aluminum for the adherends, and Metlbond 1113 and 
AF-126-2 for the adhesives. 

5.1 Materials Properties 

The mechanical and thermal properties for the materials used in 
this study were taken from the literature whenever possible. However, 
as was stated in previous chapters, complete properties were not always 
available for a given material system. Therefore, laminate properties, 
with the exclusion of uni-directional laminates, were predicted in 
accordance with the analysis procedures presented earlier and the 

results csn be found in Appendix C which contains all of the material 
properties of this study. 

5.2 Averaging of Finite Element Results 

The finite element analysis presented in Chapter 3 is based upon a 
displacement formulation. This approach yields results in the form of 
displacements at the node points and stresses and strains which are 
constant over each element. Because the stresses are constant for an 
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individual element, a distribution of stresses over a series of elements 
may appear to be discontinuous. In many cases, however, it is known 
that the stresses must be continuous. For this reason, stress averaging 
was used to produce the desired smooth distributions. 

Noting Fig. 6a, stresses are presented along the line A-A which 
corresponds to the mid-plane of the adhesive layer. The stresses pre- 
sented at point E would correspond to an average of the stresses in 

elements 3 and 4. This method of averaging was used for all joints with 
adhesive layers. 

When considering bonded joints without adhesive layers, the stresses 
in question are along the interface between the two adherends. This 
interface is shown as line C-C in Fig. 6b. For the stresses t ^ and 
which must be continuous across the interface, the results presented 
correspond to an average of elemental stresses above and below the inter- 
face. Thus, these stresses at point G would consist of an average of 
elements 13, 14, 21, and 22. The normal stress components and are 
not necessarily continuous across this interface so these stress 
components are averaged along both line B-B, and line D-D. At point F, 

an average of element 11 and 12 is presented and at point H, the stresses 
are averaged between element 19 and 20. 

The stress components and have not been mentioned as they 
do not occur either at the mid-plane of the adhesive layer in adhesive 
bonded joints, nor at the interface between the adherends in non- 
adhesive bonded Joints since the adherends in this study, when composite 
laminates, are considered to be homogeneous, orthotropic materials. If 






















these components were present it could easily be shown that t must be 

xz 

continuous across an adherend-adherend interface while would need not 
be. 

5.3 Finite Element Representations 

For the analysis of bonded joints, two finite element models were 
used. Fig's. 7 and 8 show partial plots of the finite element models of 
joints with, and without adhesive layers, respectively. Both of these 
models were generated using a mesh generator described by Bergner, Davis, 
and Herakovich [18]. in both figures the scaling of the model for the 
figure is not uniform. In Fig. 7, the aspect ratio of an element in the 
adhesive layer ranges from 2.5 at the free edge to 15 at the center of 
the adhesive layer. The aspect ratio's of the elements in the adherends 

range from 1.1 to 9.0. For the joints without adhesives (Fig. 8) all 
aspect ratios are 1.0. 

5.4 Stress Free Temperature 

Bonded joints are, in general, cured with a combination of elevated 
temperature and pressure. The maximum temperature involved in this 
process IS known as the cure temperature. The temperature at which 
curing stresses begin to form is the stress free temperature and, in 
general, the cure and stress free temperatures are not the same. The 
stress free temperature of the adhesives used in this study was chosen 
to be 270°F. This value was selected because both adhesives are epoxy 

based and 270°F was the value used in [1] for epoxy matrix material 
systems. 
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For bonded joints without adhesives, the stress free temperature 
was chosen to be 350°F. This is the value reported in [19] as the 
stress free temperature of graphite-polyimide laminates. This is ap- 
propriate as polyimides were the only composite lamiantes used in joints 
without adhesives. 

5.5 Linear Elastic Results 

This section contains the linear thermoelastic results for various 
joints. The dimensions for adhesive bonded joints are shown in Fig. 1. 
For joints without adhesives the dimensions are identical to those of 
joints with adhesives except the adhesive layer is removed. Most of the 
curves in this section were drawn by the VPI & SU computer plotter. 

In the figures that follow the superscripts M and T are used to 
differentiate stresses. Mechanically induced stresses are indicated by 
the superscript M while thermal, or curing, stresses are denoted by the 
superscript T. These are also used In combination indicating a super- 
position of mechanical and curing stresses. In some instances a curing 
stress is referred to, while the corresponding figure presents only the 
mechanical, and combined mechanical and curing stresses. The magnitude 
of the curing component can, of course, be determined by taking the 
difference of the combined, and mechanical stresses. 

5.5.1 Single Lap Joints with Adhesives 

5. 5. 1.1 [0] Graphite-Polyimide Adherends 

The adhesive stresses of a single lap joint with [0] graphite- 
polyimide adherends and Metlbond 1113 adhesive are shown in Fig. 9. The 



^/L 


Figure 9. Elastic Mechanical and Cui 
of a Single Lap Joint witi 
and Met! bond 1113 Adhesivi 


ng Adhesive Stresses 
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loading consists of a thermal increment of -200°F and an applied dis- 
placement. The resulting stresses are shown to £/L = 0.5 as the stresses 
are symmetric about this line. 

Upon examination of Fig, 9 it is seen that the stress free boundary 
conditions 


i/l = 0 


= 0 and T 


yz 


£/L - 0 


= 0 


(5.1) 


are not satisfied by the finite element solution. This is due to the 
nature of the constant stress finite elements used and the limitations 
on the maximum number of element available. In order to check the 
finite element analysis' ability to meet such stress free boundary 
conditions, an analysis was performed on a small portion of the adhesive 
layer from £/L = 0.0 to £/L = 0.05. The displacements predicted along 
the upper and lower interfaces of the adhesive in the joint solution 
were used as loading for the partial adhesive analysis. The stress 
distributions produced by this analysis exhibited the proper trends 
with 0 ^ reaching a peak value near £/L = 0 and 0 and t tending 
towards zero. These distributions are not presented, however, as they 
appeared very erratic. It is believed that this was caused by round- 
off error in the applied displacements. This error may have become 
significant after subtracting rigid body motion from the joint analysis 
displacements. 

In order to check the validity of the finite element solution 
presented in Fig. 9, a number of static equilibrium calculations were 


B 


(b) 


ee Body Diagrams of Si 
der Two Sets of Bounda 
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made. The equi nbrlum equations for one half of a single lap joint 
corresponding to the free body diagram of Fig. 10a are 


SFy = 0 = - R^, 

= 0 = 

and (5.2) 

sMq = 0 = 1 tV^ - BN^ 

The finite element program NONCOMl does not back substitute the nodal 
displacements to solve for the nodal forces. Therefore it is not 
possible to determine the reactions or However, the reaction Rj^ 
can be determined as the average of the a stresses of the elements 
adjacent to the edge where acts multiplied by the adherend thickness 
and assuming a unit depth, A comparison of V-j and determined from 
the finite element solution indicates a four percent error as shown in 
Table 2. 

Since the unknown reactions severely limit the equilibrium calcula- 
tions for the previous joint, similar calculations were also performed 
on a more simply constrained joint. This joint corresponds to the joint 
shown in Table 1, condition 4, and a free body diagram of one half of 
the joint is shown In Fig. 10b. For this free body diagram the equili- 
brium equations are 



TABLE 2 


Results of Equilibrium Checks of Single Lap 
Joint with [0] Gr/Pi Adherends and Metlbond 1113 Adhesive 
Under Two Sets of Boundary Conditions 


Joint F.B.D. 


{Fig. 10) 

y 


V-, = 92.6 lb 

a 

Rj^ = 89.0 lb 


ERROR =4.0% 


V2 = 71.8 lb 

b 

R^ = 78.3 lb 


ERROR =8.3% 


N2^ = -0.25 lb 
NEXACT ^ hj 
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/ 

and 

eHq, = 0 = ^ tVg 

As in the previous joint, the reaction moment cannot be determined 

from the finite element solution. However, by removing the constraint in 

the 2 direction, the integral of the stresses, Ng, must now equal zero. 

The results of these equilibrium calculations are also shown in Table 2. 

Returning to Fig. 9, it can be seen that under mechanical loading 

only, the shearing stresses dominate the stress fields. This is a 

function of the overlap lenth, L. If the overlap were longer, for a 

given loading, the shearing stresses would be reduced while the peak 
M 

stresses would increase. This will be shown later in the section 
containing nonlinear results and can be verified by considering the 
force and moment equilibrium equations (5.2), It is interesting to note 
that the normal stresses and are very close in magnitude for 

this joint. 

Since the stresses presented are produced by the displacements of 
the joint, its deflected shape would provide significant insight into 
the physics of the problem. Fig. 11 shows the deflection of the upper 
edge of the lower adherend under mechanical loading only. The dashed 
line signifies the beginning of the overlap {£ = 0). In this figure 
it is difficult to distinguish any curvature of the adherend in the 
region of the overlap because of the relatively small distance involved. 

A plot of the overlap only showed a nearly straight line distribution 
also indicating very little bending in this region. It is interesting 
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Figure 11. Deflected Shape of the Lower Adherend of a Single Lap Join Under 
Mechanical Loading 
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that the maximum bending occurs just before the overlap region. However, 
due to the layout of the finite element model, which was designed for 
adhesive studies, the stresses in this high bending region cannot be 
accurately obtained. 

Again returning to Fig. 9 it 1s seen that the only nonzero components 
of curing stresses are and o-^. Actually, the finite element solution 
did predict other components of curing stresses but their magnitudes 
were Insignificant. It is difficult to make a valid comparison of the 
relative magnitudes of the mechanical and curing stress components as 
the mechanical loads were produced by a small load increment while the 
curing stresses are due to the full temperature change from the stress 
free temperature to room temperature. The curing stresses represent 
approximately 15 percent of the ultimate strength of the adhesive and 
while this magnitude is not exceedingly large, its contribution should 
be Included 1n a failure theroy. 

5. 5. 1.2 [0/±45/90] Graphite-Polyimide Adherends 

F1g. 12 represents the mechanical and curing stresses for a single 
lap joint with [0/+45/90] graphlte-polyimlde adherends and Metlbond 
1113 adhesive. The loadings are identical with those of the previous 
joint. Comparing Fig's. 9 and 12 it can be seen that under mechanical 
loading only, the magnitude of the peak adhesive stresses decrease with 
decreasing adherend stiffness E„. Thus, for the same loading, the joint 

O' 

with [0] adherends has higher stresses than the joint with [0/±45/90] 
adherends. However, when considering curing stresses, this is no longer 
the case. It can be seen that the magnitudes of the stresses al and 

A 



ig Adhesive Stresses 


;0/±45/90] Gr/Pi 
Adhesive ^ 






are larger for the [0/±45/90]^ joint than for the [0] joint. This 
is due to the effects of the transverse stiffness and coefficient of the 
thermal expansion. For the [0/+45/90] adherends the longitudinal (y) 
and transverse (x) directions have identical modulus and coefficient of 
expansion. This quasi-isotropy leads to curing stresses, and 

A 

that are identical for the joint with [0/±45/90]g adherends. This 
is not the case for the joint with [0] adherends as the transverse 
direction has a much lower modulus and higher coefficient of expansion 
than the longitudinal direction. 

Because the adherends of this joint have identical transverse and 
longitudinal stiffness it is appropriate to determine what effects the 
transverse stiffness have upon the mechanically induced stresses. For 
this purpose the joint of Fig. 12 was also loaded under a classical 
plane strain assumption. The resulting stresses are shown in Fig. 13, 
The loading and materials were identical for both joints except that for 
the joint of Fig. 13 the average normal stress acting perpendicular to 
the plane was not specified and the strain normal to the plane, e , was 

A 

zero. Also, stresses are not presented for the classic plane strain 
solution. Comparing the two figures (12 and 13) it can be seen that the 
mechanically induced stresses are reduced slightly for the quasi three- 
dimensional analysis. This indicates that a 2-D solution would under- 
estimate the strength of this joint. Comparisons of the stress com- 
ponents for the two joints can be found in Table 3, The curing stresses 
are identical for the two joints because both analyses were performed 
under the 3-D analysis. 


Figure 13. Elastic Mechanical and Curing Adhesive Stresses 
of a Single Lap Joint with [0/±45/90] Gr/Pi 
Adherends and Metlbond 1113 Adhesive,^e =0 
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5. 5. 1.3 [90] Oraphite-polyimide Adherends 

The mechanical and curing stresses of a single lap joint with [90] 
graphite-polyimide adherends and Metlbond 1113 adhesive are shown in 
Fig. 14. Comparing Figs. 9, 12, and 14 the trends pointed to earlier 
concerning adhesive stresses and adherend stiffness are again confirmed. 
It is seen that the [90] adherends produce the lowest mechanical stress- 
es of the three joints. Comparing Fig's. 9 and 14 only it is revealed 
that the curing stresses for the [0] joint are identical to the o[ 

stresses for the [90] joint. This is as would be expected and the same 

T T 

correspondence is also seen between the a of the [90] joint and the a 

X y 

of the [0] joint. 

5. 5. 1.4 [0] and [90] Graphite-PolyimIde Adherends 

The adhesive stresses for a single lap joint with [0] and [90] 
graphite-polyimide adherends and Metlbond 1113 adhesive due to mechani- 
cal loading only are presented in Fig. 15. The most striking aspect of 
these stresses is the lack of symmetry present in their distributions. 
This is due, of course, to the unsymmetric nature of this joint. It is 
also interesting to note that the peak values of stress occur near the 
line z/l = 1.0 This correspondence of the peak stresses with the more 
flexible adherend seems inconsistent as the trends of the previous 
joints pointed to higher stresses with stiffer adherends. This can be 
explained upon examination of Fig. 16 which presents displacements of 
the upper adherend-adhesive interface relative to the displacements of 
the lower adherend-adhesive interface. The displacements are normalized 
with respect to the thickness of the adhesive layer. It can be seen 
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Elastic Mechanical Adhesive Stresses of a 
Single Lap Joint with [0] and [90] Gr/Pi 
flHhprpnds and Metlbond 1113 Adhesive 







that the relative displacements (y-component) reach a much higher 
value at Ji/L = 1 than at Ji/L = 0, while the displacements (z-com- 
ponent) are nearly symmetric about the center of the joint {n/L = 0.5). 
The normalized component v^ can be considered as a shear strain thus in- 
dicating why the shearing stresses are much larger near i/l = 1 in Fig. 
15. It is interesting that the component is nearly symmetric. This 
indicates that bending is not a major factor. This corresponds to the 
results seen in Fig. 11 where the curvature of adherend was nearly zero 
in the overlap region. From Fig. 16 it is seen that the increased 
stresses are due, almost entirely, to the increased flexibility (E ) of 
the [90] adherend. 

M 

Returning to Fig. 15 a discontinuity in the cr stress distribution 
can be observed at z/l = 0.2. This can be explained as a change in 
adherend stiffness corresponding to a change in the finite element 
representation of the adherend. Referring back to Fig. 7, this change 
in representation is seen as the point at which the adherend in the 
model is changed from two, to one layer of elements. This situation was 
unavoidable due to the large number of elements required to model the 
adhesive layer and a limitation on the maximum number of elements 
available. This change in stiffness Is recognizable in many of the 
stress distributions presented for adhesive bonded joints. 

Curing stresses for the single lap joint with [0] and [90] ad- 
herends are presented in Fig. 17. These stresses are, not surprisingly, 

different than those for the joints presented earlier. For this case, 
all of the components of stress induced by mechanical loading are 
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present due to cure. From observing the and distributions it can 
be seen that the integrals of these stresses appear to be zero, as would 
be required by equilibrium considerations. It is interesting to note 
the magnitude of the stress peaks near a/L = 0 and £/L = 1 . The 
values of shearing cure stresses are nearly 75 percent of the ultimate 
shear strength of the adhesive. Because the mechanical and curing 
stresses have opposite signs near i/l = 1 it can be seen that curing 
counteracts the large mechanical stresses in this region. It should be 
mentioned, however, that the percent modulus retentions, used for de- 
termining moduli as functions of temperature, are only strictly valid 
through a limited portion of the stress-strain curve of a material. 
Because of this, the peak shear cure stresses may not be quite as ac- 
curate as the other components of cure stress which correspond to points 
lower on their respective stress-strain curves. 

5.5.2 Double Lap Joints with Adhesives 

5.5.2. 1 [0] Graphite-Polyimide Adherends 

The adhesive stresses due to curing and mechanical loading for a 
double lap joint with [0] graphite-polyimide adherends and Metlbond 1113 
adhesive are presented in Fig. 18. Here, as with the single lap joints, 
the mechanically induced shearing stresses dominate the mechanical 
stresses. It is interesting to note the lack of syimietry present in 
these stress distributions due to the restrictions upon the w displace- 
ments along the midsurface of the inner adherend. These restrictions 
are induced by the symmetry of joint about the midsurface of the inner 
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Figure 18. 

Elastic Mechanical and Curing Adhesive Stresse 
of a Double Lap Joint vdth [0] Gr/Pi Adherends 
and Metlbond 1113 Adhesive 
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adherend. Comparing Fig's. 9 and 18 it can be seen that the mechanically 

induced stresses are smaller for the double lap joint than for the 

single lap joint. The distributions are more uniform for the double lap 

joint indicating reduced bending in the adherends. This is also due to 

the restrictions on the w displacements. 

The curing stresses for this double lap joint follow the trends 

exhibited by the single lap joints where the adherends were identical, 

T T 

with and being the only significant curing stresses. As expected, 
the magnitudes of the curing stress components are identical to those of 
the single lap joint with [0] adherends. 

5. 5. 2. 2 [90] Graph! te-Polyimide Adherends 

Fig. 19 represents the curing and mechanical loading induced stress- 
es of a double lap joint with [90] adherends. Comparing Fig's. 18 and 
19 it is seen that the peak stress values decrease with increased flexi- 
bility of the adherends, as was the case with single lap joints. Noting 
the curing stresses of both of these joints it can be seen that the 
cure stresses are higher for the [90] adherend joint while the aj 
cure stresses are higher for the [0] adherend joint. This is exactly 

the same as with the single lap joints and the reason for it is also the 
same. 

5. 5. 2. 3 [0] and [90] Graphite-Polyimide Adherends 
Mechanically induced adhesive stresses of a double lap joint with 

[0] and [90] graphite-polyimide adherends and Metlbond 1113 adhesive are 
represented in Fig. 20. As was the case with the two previous double 
lap joints, comparisons made with a single lap of the same adherends 
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Figure 19. Elastic Mechanical and Curing Adhesive Stresses 
of a Double Lap Joint with [90] Gr/Pi Adherends 
and Metlbond 1113 Adhesive 
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Figure 20. Elastic Mechanical Adhesive Stresses of a Double 
Lap Joint with [0] and [90] Gr/Pi Adherends and 
Metlbond 1113 Adhesive 


(MPa) 



(Fig. 15) shows lower peak stress values and a more uniform distribution 
of stresses for the double lap joint. A striking similarity between 
these two joints is the location of the peak stress values corresponding 
to a line near ifl = 1.0. In the double lap joint, as in the single 
lap, this is due to the increased relative displacement v^ at the line 
£/L = 1.0. It is interesting to note that the restrictions on the 
displacements in the z direction in the double lap do not greatly de- 
crease the relative difference in peak stresses near £/L = 0 and £/L = 
1.0 in comparison to the single lap joint, indicating that the effects 
of bending are also small for the double lap joint. 

The curing stresses for this double lap joint are presented in Fig. 
21. This figure shows that the curing stresses are very small while 
the stresses are approximately 90 percent of the ultimate shear 
strength. Comparing these stresses to those for a single lap joint of 
the same adherends (Fig. 17) indicates that for curing stresses, the 
restrictions upon displacements at the midsurface of the inner adherend 
in the double lap joint are not necessarily helpful. The double lap 
joint produces higher shearing stresses, but lower peel stresses than 
those of the single lap joint. The only difference between the two 
joints is the increased bending stiffness of the inner adherend in the 
double lap. This indicates that bending stiffness is an important 
factor in these curing stresses. This must also be the cause for 

differences in the a and a curing stresses of the two joints. 

X y 
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Before presenting the first set of results in this section it should 
be pointed out once more exactly how the following stresses were 
averaged (see section 5.2). For the case of stresses which must be 
continuous across the interface between the adherends, the stresses pre- 
sented are an average of the elemental stresses above and below this 
interface. For stresses which need not be continuous, the results 
presented consist of elemental stresses averaged either above or below 
the interface, thus two and curves are presented for each joint. 
This will also be the case for the stress results in the next section. 

5. 5. 3.1 [0] Graphite-Polyimide Adherends 

Fig. 22 presents the interfacial stresses of a single lap joint 
with [0] graphite-polyimide adherends and no adhesive. These stress 
distributions are interesting in that the shear stress is not the 
dominant stress as it was in the joints with adhesives. For this joint 
it is seen that the oj upper and lower stresses are of greater magni- 
tude and that a” attains higher peak values than the shearing stresses. 
For this joint xj^ appears to be nearly uniform along the interface at 
a relatively low value. Due to the strength of the [0] laminates in 

the fiber direction it would appear that failure would initiate as a 

M 

result of the peel stresses 

Curing stresses are not presented for this joint as they do not 
exist. The finite element solution was also checked on this point w^th 

the results being zero as required. 

5. 5. 3.2 [90] Graphite-polyimide Adherends 

Mechanically induced interfacial stresses of a single lap joint 
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with [90] graphite-polyimide adherends are presented in Fig. 23. These 
distributions show the same trends as the previous joint with the 
upper and lower stresses being the largest, reaching a relatively 
high peak, and snowing a nearly smooth distribution. Comparing 
Fig's. 22 and 23 it is seen that the larger stress values occur with the 
stiffer adherends as was the case for all of the adhesive bonded joints 
in previous sections. 

5. 5. 3. 3 [0] and [90] Graphite-Polyimide Adherends 

Fig. 24 is a plot of interfacial stresses for a joint with [0] and 

[90] graphite-polyimide adherends under mechanical loading. As with the 

previous two joints, the internal stresses upper and lower are the 

largest in magnitude. It is interesting to note that for this joint the 

peak stresses do not necessarily occur with the more flexible adherend, 

as was true for adhesive bonded joints. For this joint and 

z y 

lower have peak values near l/L - 0 while and upper have peaks 
near z/l = 1.0. Further examination of Fig. 24 indicates that the 
internal stresses are highest in the direction of the fiber in these 
unidirectional laminates. Thus is larger for the [0] adherend than 
in the [90] adherend while is larger for the [90] adherend. 

Curing stresses for this joint are presented in Fig. 25. The 
magnitudes of the stresses in the [0] adherend are nearly 60 percent 
of the ultimate strength of laminate which represents a significant 
curing stress. The a curing stresses in the [90] adherend are also 

•j 

large at 40 percent of ultimate. The other internal components of 
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Figure 23. Elastic Mechanical Interfacial Stresses of 
a Single Lap Joint with [90] Gr/Pi Adherends 
and no Adhesive 
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Figure 24. Elastic Mechanical Interfacial Stresses of 
a Single Lap Joint with [0] and [90] Gr/Pi 
Adherends and no Adhesive 
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Figure 25. Elastic Curing Interfacial Stresses of a Single 
Lap Joint with [0] and [90] Gr/Pi Adherends 
and no Adhesive 
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T T 

stress upper and lower are of comparable numerical value, but are 

in the fiber direction of their respective adherends and therefore do 

not represent such significant percentages of the ultimate strengths. 

This figure shows that equilibrium appears to be satisfied by the 
T 

and t' curing stresses. 

jf ^ 

5.5.4 Double Lap Joints without Adhesives 

Fig's. 26 and 27 represent mechanically induced interfacial stress- 
es for double lap joints with [0] and [90] graphite-polyimide adherends 
respectively. These joints show higher stresses with the stiffen adher- 
ends as with all other joints presented. It can be seen that with these 

M 

joints the a ' stresses are the largest as has been shown for all joints 
%/ 

without adhesives. 

Comparing Fig's 22 and 26 it can be seen that the double lap joint 
has higher stresses for the, same displacement loading. This can only be 
caused by the increased bending stiffness of the inner adherend in the 
double lap joint. This trend is also present in comparison of Fig's. 23 
and 27 representing joints with [90] adherends. For joints with ad- 
hesive layers the effect of the increased bending stiffness of the inner 
adherend was to smooth and reduce slightly the stress distributions for 
double lap joints in comparison to single lap joints. 

Fig's. 28 and 29 represent the mechanical and curing stresses 
respectively for a double lap joint with [0] and [90] graphite-polyimide 
adherends. As with the single lap with [0] and [90] adherends and no 
adhesive, the peak stresses for mechanical loading do not occur with the 
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Figure 27. Elastic Mechanical Interfacial 
Double Lap Joint with [90] 6r/ 
no Adhesive 
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more flexible adherend. This is again contrary to the results presented 
for lap joints with adhesive layers. The curing stress components 
for the [0] adherend and for the [90] adherend represent approxi- 
mately 60 percent of ultimate strength for the laminates. 

5.5.5 Elastic Loading Comparisons 

After reviewing the stress distributions of the four previous sub- 
sections, interesting comparisons can be made by considering the force 
loading corresponding to the displacement applied for each of the joints. 
The results can be seen in Table 4. In this table the force load for 
double lap joints corresponds to the total load carried by the joint. 

The forces for all of the joints are calculated by averaging stresses in 
the elements at the ends of the adherends, multiplying by the thickness 
of the adherend, and assuming a unit depth. 

Making comparisons of single and double lap joints with the same 
adherends it can be seen that the double lap joints carry more than 
twice as much force as the single lap joints. In the cases where the 
joints have adhesives, this can be seen as a beneficial effect of the 
increased bending stiffness of the inner adherend, as the peak stresses 
for double lap joints of this category are lower than for the single lap 
joints. The increased load carrying capacity is due to the more uniform 
shear stress distribution of the double lap joints, which creates a 
larger resultant force opposing the load. 

Comparisons of single and double lap joints without adhesives also 
show a more than doubled load carrying capacity for the double lap 
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TABLE 4 

Force-Displacement Result Comparisons for 
Elastic Joint Solutions 


iBIBM 

Adherends 

Adhesive 

Bonded 

Applied 

Displacement 

Corresponding 

Force 

Single 

[0], [0] 

Yes 

0.0005 in. 

88.8 lb. 

Single 

[90], [90] 

Yes 

0.0005 in. 

11.6 lb. 

Single 

[0/±45/90]^, 

[0/±45/90]^ 

Yes 

0.0005 in. 

38.9 lb. 

Single 

[0], [90] 

Yes 

0.0005 in. 

20.9 lb. 

Double 

[0], [0], [0] 

Yes 

0.0005 in. 

202.0 lb. 

Double 

[90], [90], [90] 

Yes 

0.0005 in. 

24.4 lb. 

Double 

[0], [90], [0] 

Yes 

0.0005 in. 

44.1 lb. 

Single 

[0], [0] 

No 

0.0005 in. 

86.3 lb. 

Single 

[90], [90] 

No 

0,0005 in. 

10.9 lb. 

Single 

[0], [90] 

No 

0.0005 in. 

20.2 lb. 

Double 

[0], [0], [0] 

No 

0.0005 in. 

244.2 lb. 

Double 

[90], [90], [90] 

No 

0.0005 in. 

34.7 lb. 

Double 

[0], [90], [0] 

No 

0.005 in. 

159.0 lb. 
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joints. However, in this case, the peak stresses are higher for the 

double lap joints making it unclear as to which joints are the more 
efficient. 

5.6 Nonlinear Results 

Bonded joint stress distributions presented in this section were 
predicted using the analysis of Chapter 3 and the nonlinear formulation 
of Chapter 4. Mechanical and thennal loading was applied as a series of 
increments. Where curing stresses are presented, they correspond to 
the total curing load. Mechanically induced stresses are presented at 
three load levels for each individual joint. The load levels for an 
individual joint do not necessarily correspond to those of other joints. 
The dimensions of the joints of this section are identical to those of 
the elastic results (Fig. 1) except where otherwise stated. 

5.6.1 Single Lap Joints 

5. 6. 1,1 Lap Shear Test 

The adhesive shear stress and strain distributions for a lap shear 
test joint with aluminum adherends and Metlbond 1113 adhesive are shown 
in Fig's. 30 and 31 respectively. This joint corresponds to that used 
In Ref. [20] for determining the adhesive shear properties as used in 
this study. The dimensions of this joint are presented in Table 5. 

These two figures (30 and 31) point to the effects of the nonlinear 
shear behavior of the adhesive. Examination of Fig. 30 reveals that as 
the displacement loading increases, the shear stress distribution 
becomes more uniform. Fig. 31, however, shows that the shear strain 
















Figure 30. 
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Lap Shear Joint wi 
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distribution does not become more ur^fiorm as loading increases. This 
can be explained by considering the shear stress-strain response of the 
adhesive (Appendix C, Fig. C.6). The slope of this curve becomes much 
smaller as the strain is increased. Therefores in the bonded joints the 
increment of stress corresponding to an increment of strain, at high 
strain lev ’’s, is smaller than at low strain levels. Thus, the shear 
stress distribution in the center of the joint is increasing more 
rapidly as loading increases than at the edges {z/L = 0 and = 1 ) , 
because the strains are higher at the edges. It is interesting to note 
that the largest shear stress presented (Fig. 30) is nearly a constant 
value. This value corresponds to the ultimate shear strength of the 
adhesive. 

As was stated earlier, this joint corresponds to a joint used in 
the lap shear tests [20]. Therefore, comparisons between numerical and 
experimental work were made and the results can be found in Table 6 and 
Fig. 32. 

The experimental stresses and strains presented in Table 6 are 
values corresponding the maximum stress of the adhesive shear curve. 
These stress and strain values are chosen for the comparison because of 
the nature of both the Ramberg-Osgood [17] approximations, and the 
finite element analysis. The Ramberg-Osgood parameters cannot model the 
stress-strain curve beyond the point at which the slope becomes zero and 
the finite element formulation does not produce a positive definite 
stiffness matrix when a negative modulus is used. 

In Table 6, the numerical strain chosen for the comparison was near 


Max Adhesive 
Shear Strain 
{%) 


Max Adhesive 
Shear Stress 
{ks1) 


Failure Load 
(kips) 


Numerical 27 


Experimental 29* 

[ 20 ] 


Percent 7 % 

Difference 



Corresponds to maximum stress value 




1.35 

1.36 

0.7% 


















the line £/L = 0. Thir> is also the region where the strain was de- 
termined in Ref, [20], This table shows good correlation between ex- 
perimental results and the numerical prediction. In Fig. 32 two finite 
element stress-strain responses are presented. In both of these finite 
element curves ^ the stress corresponds to an average of elemental 
stresses throughout the adhesive layer. 

The upper numerical curve (AVE.F.E.) presents strains that are also 
averaged over the entire adhesive layer. The lower numerical curve 
{F.E. STRAIN NEAR 2 ,/L = 0, AVERAGED STRESS) shows strains corresponding 
to the finite elements adjacent to the line z/l - 0 in the adhesive. 

The upper curve shows a better stress correlation while the lower curve 
shows a better ultimate strain correspondence with experimental data. 

Fig's. 33, 34, and 35 present the aj), and adhesive stresses 
of the same lap shear test joint under identical loadings. These 
figures do not show the effects of adhesive nonlinearity in such a 
pronounced fashion as Fig, 30, The reason for this is two fold. First, 
the extensional stress-strain response of the adhesive is not as non- 
linear as the shear response. Secondly, the extensional stress values 
produced are not as large in magnitude relative to the ultimate strength. 
The maximum o and o stresses correspond to approximately 15 percent of 
the ultimate while the maximum a component is approximately 10 percent 
of ultimate. Upon examination. Fig's. 34 and 35 reveal the discon- 
tinuity at l/^ = 0.2 that v/as discussed in the elastic results (section 
5. 5. 1.4). Thermal stresses are not presented for this joint as they 
would have no bearing on the comparisons made. This is because the 


Figure 33. Nonlinear Mechanical Adhesive Stresses of a 
Lap Shear Joint with Aluminum Adherends and 
Metlbond 1113 Adhesive 






Figure 34. Nonlinear Mechanical Adhesive Stresses of a 
Lap Shear Joint with Aluminum Adherends and 
Met! bond 1113 Adhesive 
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Lap Shear Joint with Aluminum Adherends and 
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shear properties produced from the experimental work on this joint [20] 
included any effects of curing. 

An interesting comparison can be made concerning results obtained 
from the two different solution formulations presented in Chapter 3. In 
order to produce similar stress and strain distributions, it was neces- 
sary to apply a larger displacement load to the joint analyzed under the 
quasi 3-dimensional formulation. The additional load corresponded to 
1/15 of the total displacement applied for the 2-dimensional solution. 
The resulting ultimate force loads were identical for the two analyses 
indicating that under the quasi 3-dimensional formulation the joint was 
more flexible and therefore capable of withstanding a larger displace- 
ment load. 

5,6, 1.2 [0/90/0/90/0] Boron-Epoxy Adherends and 
AF-126-2 Adhesive 

This adhesive bonded joint was selected for analysis in order to 
compare the results of this study to those of Ref. [11]. For this 
comparison it was necessary to use a force loading instead of the dis- 
placement loading used for all other joints in this study. The force 
loading was required in order to exactly match the loading in [11]. In 
[11], two solution procedures are used. In the first, an iterative 
finite element analysis is used to account for material nonlinearities. 
The second procedure utilizes direct numerical integration of the 
governing differential equations for the joint. Results obtained by the 
second procedure are labeled as theoretical in the following figures. 

The dimensions of this joint are given in Table 5. 
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F1g s. 36 and 37 present the adhesive stress components and 

M . 

respectively. In Fig. 36 it can be seen that the distributions of 
the shear stresses predicted by the present study compare favorably with 
the two distributions predicted in [11]. The principal difference is 
seen to be the magnitude of the peak stresses near the lines £/L = 0 and 
VL - 1. It is interesting that none of the shear stress distri- 
butions presented in Fig. 36 satisfy the stress free boundary condition 
^yz " ^ • It appears that the areas under each of the three 

curves are approximately the same as required by equilibrium considera- 
tions. 

In Fig. 37 it can be seen that the distributions predicted by 
the present study and theoretical results of [11] have significant 
differences. The present analysis predicts a symmetric distribution of 
peel stresses while the theoretical results [11] show peak stresses at 
A/L = 0 and i/l = 1 differing by more than 100 percent. The nonsym- 
metric nature of these results appears physically inconsistent as the 
adherends are identical. The finite element results presented in [11] 
do show the symmetry of stresses as predicted by this study. Another 
interesting aspect of the theoretical results is the reversal in sign of 
the peel stresses near i/l = 1.0. This is not seen in the finite 
element results of this study or [11]. 

Comparisons of the numerical values of the peel stresses cannot 
realistically be made as the extensional properties of the AF-126-2 
adhesive were not known. The extensional Ramberg-Osgood coefficients 
used for this adhesive correspond to the extensional properties of 
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Figure 36. Nonlinear Mechanical Adhesive t Stresses of a 
Single Lap Joint with [0/90/0/90/0] B/E 
Adherends and AF-126-2 Adhesive 
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Figure 37. Nonlinear Mechanical Adhesive Stresses of a 
Single Lap Joint with [0/90/0/90/0] B/E 
Adherends and AF-126-2 Adhesive 
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Metlbond 1113. 

Comparing Fig's. 36 and 37, it can be seen that for all solutions 
except the finite element solution of [11], the peel stresses are 
larger than the shear stresses. This verifies the statement made about 
the effects of the overlap length, L, upon the magnitudes of the shear- 
ing and peel stresses in an earlier section (5. 5. 1.1). For this joint, 
with a large adherend overlap, the peel stresses are dominant with 
respect to maximum value. If the magnitudes were compared with respect 
to ultimate strengths however, it is believed that the shearing stresses 
would again dominate. It is not known what these ultimate strenths are 
however, and therefore this comparison cannot be made. 

The curing stresses for this joint are presented in Fig. 38. As 

was shown for the elastic results, the only significant curing stresses 

are cr^ and This is again due to the identical adherends and the 

material properties of the adherends and adhesive. It is interesting 

that even though the values of the curing stress components and 

2 yz 

are insignificant, they reach peak values nearly an order of magnitude 
larger than in any other joint with identical adherends. Since these 
curing components and are negligible, the stresses presented in 
the previous two figures {36 and 37) can be considered either mechanical 
or combined mechanical and curing stresses. 

The magnitudes of the curing components and aJ with respect to 
the ultimate extensional strength of the adhesive is not known because 
as was stated earlier, this ultimate strength is not known. The analy- 
ses of [11] ignore the effects of curing, and while this does not affect 
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Figure 38. Nonlinear Curing Adhesive Stresses of a Single 
Lap Joint with [0/90/0/90] B/E Adherends and 
AF-126-2 Adhesive 
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the a and t stress components, it can be seen that curing is signi- 
z yz 

ficant in the two components of stress not presented in [11] (a^ and 

Oy). 

Fig's. 39 and 40 present the combined mechanical and curing stress- 
es, and of this joint. These distributions reveal that the 
y X 

curing stresses (Fig. 38) are of the same sign as the mechanically in- 
duced stresses and therefore would be detrimental to the performance of 
the joint under tensile loading, 

5. 6. 1.3 [0] Graphite-polyimide Adherends and Metlbond 

1113 Adhesive 

Nonlinear curing stresses for this joint are not plotted as the 

significant components (a^ and o^) are uniform throughout the adhesive. 

X y 

The numerical values for these stresses can be found in Table 7. This 
table presents a comparison between elastic and nonlinear results for 
this joint. For these curing stresses, the elastic solution under- 
estimates the and aj components by 10 percent and 5 percent respec- 
tively. 

The mechanical loading of this joint was analyzed utilizing the 2- 
dimensional formulation (Chapter 3). This was done because for this 
joint, the two formulations produce negligible differences in the stress 
components a\ and . The similarity in results is due to the 
relatively small difference in magnitudes of the adherend transverse 
stiffness (E ) and the adhesive extensional stiffness. Combined mech- 
anically and curing induced adhesive stresses are presented in Fig's. 

41, 42, and 43. As in the lap shear joint (Fig. 30) only the adhesive 
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Figure 39. Nonlinear Mechanical and Curing Adhesive 

Stresses of a Single Lap Joint with [0/90/0/90/0] 
B/E Adherends and AF-126-2 Adhesive 
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TABLE 7 


Comparison of Elastic and Nonlinear Adhesive Stresses for 
Single Lap Joint with [0] Gr/Pi Adherends and Metlbond 1113 Adhesive 
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Figure 43. Nonlinear Mechanical and Curing Adhesive 

Stresses of a Single Lap Joint with [0] Gr/Pi 
Adherends and Metlbond 1113 Adhesive 
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shearing stresses for this joint (Fig. 41) exhibit nonlinear effects. 
Comparisons between elastic and nonlinear mechanically induced stresses 
are presented in Table 7. These comparisons show that the and 
stress components behave linearly throughout the range of mechanical 
loading applied. The shearing stress is decidedly nonlinear though, 
with the elastic solution predicting stresses 82 percent higher than the 
nonlinear results indicate, 

5. 6. 1.4 [±45]g Graphite-Polyimide Adherends and 
Met! bond 1113 Adhesive 

The nonlinear curing stresses for this joint are not plotted for 
the same reasons as the previous joint. The numerical values of the two 
significant curing stresses are presented in Table 8 which presents a 
comparison of elastic and nonlinear results for this joint. This 
comparison shows a six percent increase in the curing stresses for the 
nonlinear analysis. 

Combined mechanical and curing adhesive stresses for this joint are 
shown in Fig's 44, 45, 46 and 47. Again, it is seen that nonlinear 
behavior is present in only the shear stresses. It is unfortunate that 
the maximum loading for this joint did not produce a peak shear stress 
corresponding to the ultimate strength. This would have produced much 
more pronounced nonlinear effects. Comparisons between elastic and 
nonlinear mechanical adhesive stresses are also presented in Table 8. 

It is seen that at the maximum load level attained, the nonlinear 
results predict a peak shear stress 28 percent below the elastic re- 
sults. The elastic mechanical results are determined as the first 
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TABLE 8 


Comparison of Elastic and Nonlinear Adhesive Stresses for a Sinai 
Lap Joint with [±45]^ Gr/Pi Adherends and Metlbond 1113 Adhesive 


Type of 

Curing Stresses 

Peak Mechanical Stresspc 

Analysis 

SET = 270°F (ksi) 


(ksi) 



"X “y 


°y 

<^2 


Elastic 

1.5 1.5 

1.3 

2.6 

2.9 

5.1 

Nonlinear 

1.6 1.6 

1.3 

2.6 

2.9 

4.0 
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increment of the nonlinear results scaled to the maximum displacement 
load. It is somewhat surprising that the nonlinear analysis does not 
predict lower stress components a”, a”, and a” for this and the previous 
joint. It can be explained however by considering the final magnitudes 
of these stress components. All are relatively low on the extensional 
stress-strain response of the adhesive. At these stress levels the 
curve is very linear resulting in very linear stress predictions. 

5.6.2 Double Lap Joints 

5.6.2. 1 Titanium Adherends Metlbond 1113 Adhesive 

Comparisons of elastic and nonlinear curing stresses for this joint 

are shown in Table 9. Once more, the two nonzero components of curing 
T 

stress and a^) were uniform and they are not plotted. Here again, 
an increase in the nonlinear results over the elastic case is seen. 

Combined stresses for this joint are presented in Fig's. 48, 49, 50 
and 51. In these figures only the shear curves show pronounced non- 
linearities, as before, but Table 9 indicates that the other stresses 
and a^) are slightly reduced for the nonlinear analysis. The 
reason this joint should exhibit nonlinear behavior where the previous 
joints did not is not immediately discernable. The most obvious dif- 
ference between this joint and the previous ones is the restriction 

placed upon the w displacements at the midplane of the inner adherend by 
the symmetry of the double lap. 

5.6.2. 2 [0] and [90] Graphite-Polyimide Adherends and 
Metlbond 1113 Adhesive 

Curing stresses for this joint are presented in Fig. 52. These 


TABLE 9 


Comparison of Elastic and Nonlinear Adhesive Stresses for a Double 
Lap Joint with Ti Adherends and Metlbond 1113 Adhesive 


Curing Stresses 
SFT = 270°F (ksi) 
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Figure 50. Nonlinear Mechanical and Curing Adhesive a 
Stresses of a Double Lap Joint with Titanium 
Adherends and Metlbond 1113 Adhesive 















Figure 52. Nonlinear uring Adhesive Stresses of a Double 
Lap Joint with [0] and [90] Gr/Pi Adherends 
and Metlbond 1113 Adhesive 
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distributions appear very similar to the elastic curing stresses for the 

same joint (Fig. 21), Comparisons between these two analyses are shown 

1n Table 10. It can be seen that the nonlinear curing stresses differ 

from the elastic results between 2 percent and 6 percent. The 23 

percent difference shown for near a/L = 0 is probably exagerated by 

the low magnitudes of these stresses. It can be seen in Table 10 that 

the nonlinear effects cause an increase in the and curing stresses 

as shown in previous results. However, the ol and stresses are 

z yz 

decreased by the nonlinear behavior. These two stresses are related 
through equilibrium (section 5. 5. 1.1) and it is therefore appropriate 
that both change in the same fashion. This decrease is due to 
the high values of shear stress produced by curing and the low adhesive 
modulus at these values. 

Fig's. 53, 54, 55, and 56 present combined mechanical and curing 
adhesive stresses for this joint. Comparing Fig's. 52 and 53 it can be 
seen that the curing stresses are very beneficial to the performance of 
this joint. In Fig. 53, as the displacement load level increases, the 
shear stresses near z/L = 1 are seen to increase more rapidly than at 
z/l = 0. Thus, the shear stresses near z/l = 1 would be much larger 
than the shear stress near £/L = 0 for mechanical loading only. Ho\'- 
ever, the curing shear stresses near z/l = 1 have the opposite sign of 
the mechanical shear stresses in this region and are of relatively large 
magnitude. Therefore, a large portion of the mechanically induced shear 
stresses are negated by the curing shear stress and thus the joint is 
capable of carrying an increased load. 


TABLE 10 


omparison of Elastic and Nonlinear Adhesive Curing Stresses for a 
Double Lap Joint with [0] and [90] Gr/Pi Adherends 
and Metlbond 1113 Adhesive 


Type of 
Analysis 


Curing Stresses SFT = 270° (ksi) 


£/L=0 


H/L=l 




a/L=l 


«./L=0 I £/L=1 


yz 

Jl/L=0 I £/L=1 


Elastic 


1.47 


1.03 


1.26 


0.67 


0.092 


-0.53 


3,6 


-4.0 


Nonlinear 


1.53 


1.08 


1.32 


0.71 


0.12 


-0.51 


3.53 


-3.8 


IDifference 


4 % 


h% 


6 % 


23 % 


4% 


2 % 


5% 





























Figure 53. Nonlinear Mechanical and Curing Adhesive t 
Stresses of a Double Laop Joint with [0] 
1113 Adhesive 
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Figure 56. Nonlinear Mechanical and Curing Adhesive 

Stresses of a Double Laop Joint with [0] and [90] 
Gr/Pi Adherends and Metlbond 1113 Adhesive 
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Chapter 6 

SUMMARY AND CONCLUSIONS 

The present analysis has been concerned with the nonlinear analysis 
of bonded joints. Upon reviewing the results presented, the following 
conclusions can be made with respect to the joint materials and geometries 
studied. 

1. The effects of adhesive nonlinearities greatly influence 
the shear stress predictions in the adhesive layer of 
bonded joints. 

2. The effects of adhesive nonlinearities have little in- 
fluence upon the normal stress components in the 
adhesive layer of bonded joints. 

3. Adherend nonlinear behavior has little effect upon 
the adhesive stresses in bonded joints. 

4. Residual curing stresses are significant in adhesive 
bonded joints. These curing stresses are detrimental 
in joints with similar adherends, but may be beneficial 
in joints with differing adherends. 

5. Residual curing stresses are significant in bonded 
joints with differing adherends and no adhesive. 

6. Residual curing stresses are not significantly in- 
fluenced by material nonlinearities or temperature 
dependent properties. 

7. Adherend stiffness has profound effects upon mechani- 
cally induced stresses in bonded joints. Stresses 
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produced in the adhesive layer of adhesive bonded 
joints and along the adherend-adherend i?iterface in 
bonded joints without adhesives are higher for 
stiffer adherends. In adhesive bonded joints with 
differing adherends, maximum adhesive stresses 
correspond to the more flexible adherend. 

8. Adhesive and interfacial stresses are non-uniform 
with maximum values produced near the edges of the 
overlap region. 

9. Adhesive bonded double lap joints are more ef- 
ficient single lap joints due to a more uniform 
stress distribution in double lap joints. 

10. A quasi 3-dimensional analysis predicts a more 
flexible joint response than a 2-dimensional formu- 
lation. A larger displacement load is required in 
a quasi 3-dimensional analysis to predict adhesive 
stresses comparable to those of a 2-dimensional 
analysis. 

11. A quasi 3-dimensional analysis demonstrates the 
effects of adherend transverse stiffness and thermal 
coefficient of expansion upon the residual curing 
stresses. 

12. Adhesive stresses are not significantly influenced 
by different symmetric loadings. Force loadings 
produce results similar to displacement loadings. 
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13. The method of solution presented satisfies static 
equilibrium very closely. 

This analysis has shown that future areas of study might include 
the following: 

1- Analysis capability for out of plane bending and 
warpage. 

2. Better representation of stress-strain response 
as a function of temperature, and moisture. 

3. Nonlinear analysis of the effects of moisture in 
bonded joints. 

4. Inclusion of a capability to allow Poisson Ratios 

to vary as a function of strain, temperature, and 
moisture. 

5. Allowing failure strengths to vary as functions of 
temperature and moisture. 

6. Consistent modeling of the interactions of tempera- 
ture and moisture. 
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APPENDIX A 

ELEMENTAL STIFFNESS MATRIX 

Equations (A.l) represent the equilibrium equations for applied 
strain loading. Equ's (A. 2) represent the equilibrium equations in 
average force loadings. In these equations, [K] is the symmetric 
elemental stiffness matrix, and {T} are force vectors correspond- 

ing to the applied strain and temperature change respectively, IF) is 
the vector of applied forces, and {x> is the vector of unknown nodal 

displacements. 

(9x9) (9x1) (9x1) (9x1) 

(A-2) 

(10x10) (10x1) (10x1) (10x1) 

Defining the following terms 

a = 

b = (Y3-Y2)/2 
c = (Z3-Z^)/2 
d = (Y^-Y3)/2 
e = {ZyZ^)IZ 
g = (Y 2 "Y ^)/2 
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A = the area of element (£) 

* 

F = average normal force 

where Y-j through and through 2^ are the coordinates of the nodal 

points of element i in the Y-Z plane, the element of the matricies of 
Equ. (A.l) can be defined as follows. 

+ Cgga2)/A^ K22 = (E55CI2 + CggC2)/A^ 

K^2 = (Cjsbd + Cggac)/A^ Kg3 = (Cg^dg + CggCe)/A^ 

'^24 “ ^^ 26 ®® * ^^5^ti)/A‘'' 

K^4 = (Cgga2 + E^gb2)/A^ K^g = (CggC^ + C^gd2)/A^ 

K^g = (Cggca + C^gbd)/A^ Kgg = (Cggec + C^gdg)/A^ 

'^16 = ^^ 26 ®® C 4 gbg)/A^ = {Cggbc + C 4 gda)/A^ 

'^17 ° ^^ 36 '^® * C^5ba)/A^ Kgg + (Cggdc + C dc)/A^ 

K^g = (Cggda + C^gbc)/A" Kgg = (Cgggc + C4gde)/A^ 

Ki9 = (Cggga + C^gbe)/A’^ 

S3 “ S4 ” (^22®^ ’’’ 

K34 = (Cggae + C^ggb)/A*' K^g = (Cggac + C^^bdl/A^ 

Sb = (Sb®® + Sb = (^2®® 

Sb = ^Sb®^ + C^ggS/A^ K ^7 = (C^^ba + Cggab)^^ 

^■57 ” C;,r-ga)/A K-n = (C«y,bc + C„-ad)/A 





where f‘s are nodal forces. 


For Equ's. (A. 2) the previously defined terms apply plus the fol 
lowing additional terms 




*^210 '' 16 '^ 


’^310 “ '' 16 ® 


^410 ° '^12® 


*^510 ° '' 12 ® 


*^610 '^ 16 ® 


K71Q = 


'^810 ' ''13'^ 


'^giO ° ''13^ 


•^lOlO " 


^10 = 5"x Fio = F* 


Tl = 

^^16 

®x ^26®y ■*■ ^36 

^2 = 


j ? T 
'^26®y 

+ ^36®! ■ 

tt 

^^16®x 

+ ^26®J 

+ '^36®! ■ 

II 


+ CggEy 

■*■ ^23 ■ 

II 

(^12®x 

+ ^22^y 

^23®z 


T 


T 


T 


'26’^xy^® 
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Tg = (C,2^x ^22^1 + ^23^2 + '^26\y^^ 

h = + ‘^23'^J + Ss^I + 

'^8 ° ^^3^x ^23®y + ^3^ + ‘^36^^"' 

^9 = (C^3"I + <^23^ * Ss"! ■ ^36^xy)9 

^10 = + =16^5^* 

where 

ej[ = (m^a^ + n^a^)^’’’ 

^y ~ (n^a^ + 111^02) AT 

4 = “3*^ 

yj[y = 2mn(a^-a2)AT 

For moisture analysis the vector {T) is identical except o.^ , ccg and a 
are replaced by 6-,, and B3. 


T" 
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APPENDIX B 

ADDENDA TO RFERENCE [1] 


During the process of modifying the program developed by Renieri 
and Herakovich [1] it was discovered that a coding error had been made 
in ■'■he equations for the transformation matrices. (Chap. 3, Equ's. 
3.9.) This was corrected and the nonlinear stress-strain response pre- 
dictions of Ref. [1] were regenerated using the same data and finite 
element model. In the majority of cases the differences proved to be 
negligible. For the cases where the differences were significant, 
comparisons between results obtained from the corrected program and 
experiment are presented. It should be noted that the coding errors 
applied only to the nonlinear results presented in Ref. [1] and not the 
elastic results. It should also be noted that no attempt to define 
failure was made by this investigator and that the last point plotted 
does not necessarily correspond to failure. 
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appendix c 

MATERIAL PROPERTIES 


This appendix contains all of the material properties used for this 
study. The data presented represents typical data from the literature. 
References are provided where appropriate. Fig's. C.l through C.7 repre 
sent the stress-strain response of the materials used. Table C.l con- 
tains the Ramberg-Osgood coefficients for these materials and Table C.2 
contains the temperature dependent properties. 

In Table C.l the symbol a* refers to the stress at which the Ram- 
berg-Osgood coefficients and become applicable. 
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TABLE C.l Ramberg-Osgood Coefficients 


El astic 
Curve Modulus 


Elastic 

Limit 

(KSI) 


n-| Ki (PSl""!) 


ng <2 (PSr"2) 


Ultimate! 

Stress 

(KSI) 



















































TABLE C.l 


Cont. 


r— 

R3 

s: 

Curve 

Elastic 

Modulus 

(MSI) 

V 

Elastic 

Limit 

(KSI) 

"1 

K-i (psr"i) 


B 

K2 (Psr"2) 

Ultimate 

Stress 

(KSI) 

1— 1 
Q. 

CC 

CD 

to 
1 — 1 
tn 

iL 

UWM 

2.95 

^^3 

1.2 

4.064 

1.856 X 10-20 




30.1 

WM 

2.95 


1.2 

4.064 

1.856 X 10-20 




30.1 

IMl 

2.95 


1.2 

4.064 

1.856 X 10-2° 




30.1 

BW 

2.95 


1.2 

4.064 

1.856 X 10-20 




30.1 

Ezz 

Tension 

2.2 

.21 

(^xz) 







9.7 


1.8 

.21 

(^xz) 


1.253 

1.9509 X 10-0 

16.8 

2.0661 

7.181 X 10-12 

21.8 

%2 

0.85 


0.1 

3.411 

2.7474 X 10-10 

12.6 

5.3158 

4.254 X 10-24 

15.0 

Gxz 

0.85 


0.1 

3.411 

2.7474 X 10-16 

12.6 

5.3158 

4.254 X 10-24 

15.0 

^xy 

5.1 








15.0 



























































[0/±45/90]s GR/PI 

































TABLE C.l Cont. 


p*~ 

+J 

fo 

s: 

Curve 

Elastic 

Modulus 

(MSI) 

V 

Elastic 

Limit 

(KSI) 

D 

K-| (psr"i) 

a* 

(KSI) 

ri2 

Kg (psr*^2) 

Ultimate 

Stress 

(KSI) 

Unidirectional B/E* 

txx 

Tension 

29.6 


1.0 

4.463 

1.8534 X 10-27 




200.5 


29.6 


1.0 

4.463 

1.8534 X 10-27 




350.0 


2.75 


1.0 

2.541 

8.1604 X 10-14 




9.6 

^yy 

Comp. 

2.75 

.35 

(vyz) 

1.0 

2.541 

8.1604 X 10-14 




45.0 

Ez2 

Tension 

2.75 

(vxz) 

1.0 

2.541 

8.1604 X 10-14 




9.6 

Ezz 

Comp. 

2.75 

. 35 
(''xz) 

1.0 

2.541 

8.1604 X 10-14 




45.0 

Gyz 

0.933 


0.1 

2.991 

7.8797 X 10-15 




9.5 

Gxz 

0.933 


0.1 

2.991 

7.8797 X 10-15 




9.5 

Sxy 

0.933 


0.1 

2.991 

7.8797 X 10-15 




9.5 


* Ref [n]. 















Curve 

Elastic 

Modulus 

(MSI) 

^xx 

Tension 

18.9 


^yy . 
Tension 

13.5 

1^91 

13.5 

Ezz 

Tension 

2.75 

BH 

2.75 

Gyz 

0.933 

Gxz 

0.933 

^xy 

0.933 


Elastic 

Limit 

(KSI) 

"1 

K-, (psr"i) 

5.16 

3.351 

6.0822 X 10-21 

5.16 

3.351 

6.0822 X 10 ' 2 '> 

5.44 

3.05 

6.7331 X 10-19 

5.44 

3.05 

6.7731 X 10-19 

1.0 

2.45 

8.1604 X 10-14 

1.0 

2.45 

8.1604 X 10-14 

0.1 

2.991 

7.8797 X 10-15 

0.1 

2,991 

7.8797 X 10-15 

0.1 

2.991 

7.8797 X 10-15 


Kp (psr"2) 


Ultimate 

Stress 

(KSI) 


129.0 


129.0 


90.7 


90.7 


9.6 


45.0 


9.5 


9.5 


9.5 

























































Metlbond 


TABLE C.l Cont 



K2 (PSr'^2) 


Ultimate 

Stress 

(KSI) 























































r-" 

4-» 

rtJ 

2: 

Temp. 

op 

% 

Exx 

% 

E.v.y 

% 

Ezz 

Unidirectional 

GR/PI 

0.0 

100.0 

100.0 

100.0 

60.0 




400.0 




500.0 




550.0 

115.0 

50.0 

50.0 

>— i 
Q- 

ct: 

CD 

1 — J 
LO 

+1 
1 1 

0.0 



100.0 

60.0 




70.0 

100.0 

100.0 


400.0 

94.0 

94.0 


500.0 

78.0 

78.0 


550.0 



50.0 


* Ref. [19]. 


** Ref. [23]. 


TABLE C.2 Thermal Properties 


% 

Gyz 

% 

Gxz 

% 

Gxy 

a-, 

(in/in/°F) 

a2 

(in/in/°F) 

“3 , 

(in/in/°F) 




0.0* 

14.1x10"®* 

14.1x10"®* 

100.0** 

100.0** 

100.0** 




93.0** 

93.0** 

93.0** 




75.0** 

75.0** 

75.0** 







0.0* 

14.1x10"®* 

14.1x10"®* 






14.1x10"® 

100.0 

100.0 







100.0 

1.7x10"® 

1.7x10"® 


93.0 

93.0 

107.0 

1.04x10"® 

1.04x10"® 


75.0 

75.0 

in.o 

.864x10"® 

.864x10"® 







14.1x10"® 
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TABLE C.2 Cont. 


Temp. 

op 

% 

Exx 

% 

% 

Ezz 

% 

Gyz 

% 

Gxz 

% 

^xy 

a-j 

(1n/1n/°F) 

ag 

(in/in/°F) 

“3 

(in/1n/°F) 

0.0 



100.0 






14.1x10-6 

60.0 




100.0 

100.0 





70.0 

100.0 

100.0 




100.0 

1.7x10-® 

1.7x10-® 


400.0 

105.0 

105,0 


93,0 

93.0 

105.0 

1.04x10-® 

1.04x10-® 


500,0 

106.0 

106,0 


75.0 

75.0 

106.0 

.864x10-® 

.864x10-® 


550.0 



50.0 






14.1x10-6 

0.0 

103.0 

103.0 

103.0 

103.0 

103.0 

103.0 

15.5x10-® 

15.5x10-® 

15.5x10-6 

100.0 

98.0 

98.0 

98.0 

98.0 

98.0 

98.0 

16.5x10-® 

16.5x10-6 

16.5x10-6 

200.0 

90.0 

90.0 

90.0 

90.0 

90.0 

90.0 

18.0x10-6 



300.0 

68.0 

68.0 

68.0 

68.0 

68.0 

68,0 

20.25x10-® 

20.25x10-6 

20.25x10-6 

350.0 

50.0 

50.0 

50.0 

50.0 

50.0 

50.0 

22.0x10-® 

22.0x10-6 

22.0x10-6 
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TABLE C.2 Cont. 



Temp. 

% 


on 

i 

Exx 


0.0 

100.0 


100.0 

100.0 


200.0 

100.0 


2B0.0 

100.0 


300.0 

100.0 


0.0 

102.0 


100.0 

100.0 


200.0 

■ 98.0 


250.0 



300.0 

97.0 


% 

Eyy 

% 

Ezz 

% 

Gyz 

% 

<^xz 

% 

^xy 

“1 

(in/in/°F) 

“2 

(in/1n/°F) 

“3 

(in/in/°F) 

120.0 

120.0 

104.0 

104.0 

104.0 

2.2x10-® 

9.6x10“6 

9.6x10-6 

92.0 

92.0 

97.0 

97.0 

97.0 

2.32x10-® 

1.11x10-5 




83.0 

83.0 

83.0 

2.4x10-6 

1.34x10-5 

1.34x10-5 



60.0 

60.0 

60.0 




45.0 

45.0 

45.0 

45.0 

45.0 

2.7x10-6 


1.68x10-5 

103.0 

120.0 

104.0 

104.0 

104.0 

2.91x10-6 

3.45x10-6 

9.6x10-6 

100.0 

92.0 

97.0 

97.0 

97.0 

2.99x10-6 

3.51x10-6 

1.11x10-5 

96.0 


83.0 

83.0 

83.0 

3.04x10-6 

3.55x10-6 




60.0 

60.0 

60.0 




93.0 

45.0 

45.0 

45.0 

45.0 

3.26x10-6 

3.71x10-6 

1.68x10-6 
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TABLE C.2 Cont. 


4J 

(O 

S 

Temp. 

Bi 

% 

^yy 

m 

N 

N 

% 

Syz 

% 

Gxz 

% 

^xy 

“1 

(in/in/°F) 

“2 

{1n/in/°F) 

'"3 

(in/in/°F) 

Titanium* 

6A1-4V 

0.0 

101.3 

101.3 

101.3 

101.3 

101.3 

101.3 

4.7x10-6 

4.7x10-6 

4.7x10-6 

70.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 




400.0 

95.0 

95.0 

95.0 

95.0 

95.0 

95.0 


5.2x10-6 

5.2x10-6 

800.0 







5.6x10-6 

5.6x10-6 

5.6x10-6 

A1 uminum* 
2024-T4 

0.0 

104.0 

104.0 

104.0 

104.0 

104.0 

104.0 

12.1x10-6 

12.1x10-6 

12.1x10’6 

70.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 




200.0 

99.0 

99.0 

99.0 

99.0 

99.0 

99.0 

12.9x10-6 

12.9x10-6 

12.9x10-6 

250.0 







13.0x10-6 

13.0x10-6 

13.0x10-6 

300.0 

94.3 

94.3 

94.3 

94.3 

94.3 

94.3 

13.1x10-6 

13.1x10-6 

13.1x10-6 

400.0 

88.7 

88.7 

88.7 

88.7 

88.7 

88.7 

13.2x10-6 

13.2x10-6 

13.2x10-6 

500.0 

80.2 

80.2 

80.2 

80.2 

80.2 

80.2 

13.5x10-6 

13.5x10-6 

13.5x10-6 
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NONCOMl FORTRAN USERS GUIDE 

Cards 1-5 (20A4) 

Column Contents 


1-80 


Title Cards 

Card 6 

(616) 


Column 


Contents 

1-6 

NE 

= Number of elements 

7-12 

NDS 

= Number of nodes 

13-18 

MDIFM 

= Number of different materials 

19-24 

NANG 

= Number of different angles 

25-30 

lELET 

s Operating temperature indicator 



0 for 70° 



> 0 for any other temperature 

31-36 

lELEM 

= Operating moisture content indicator 



0 for 0% moisture 



> 0 for elevated moisture content 

Card 7 

(616) 


Column 


Contents 

1-6 

NLOADS 

= Number of load cases 

7-12 

NPSS(l) 

= Load type number 1 



13-18 

NPSS{2) 

= Load type number 2 


etc. 

Repeated NLOADS times 

tf 


NPSS{J) 

= 1 for axial strain 


= 2 for thermal 
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= 3 for axial force 
= 4 for hygroscopic 

Card 8 (516) 


Column 


Contents 

1-6 

NINCRT(I) 

= Number of load increments for load case 

7-12 

N1NCRT(2) 

= Number of load increments for load case 

etc. 

Repeated NLOADS times 

Card 9 

(516) 


Column 


Contents 

1-6 

KEYO) 

= Print indicator for grid 

7-12 

KEY{2) 

= Print indicator for strains 

13-18 

KEY(3) 

= Print indicator for stresses 

19-24 

KEY (4) 

= Print indicator for equivalent stresses 

25-30 

KEY(5) 

= Print indicator for displacements 
KEY(I) = 0 for printing 

Card 10 

, (1016) 

1=1, NLOADS 

Column 


Contents 

1-6 

LlNCPR(l.I) 

= First increment of load case I to print 
stresses, strains, and displacements 

7-12 

L1NCPR(2,I) 

= Last increment of load case I to print 

etc. 

Repeated NLOADS times 

Card 11 

, (2F12.6) 


Column 


Contents 

1-12 

SMY 

= Scale factor for Y-coordinates 


^ I . lu 
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13-24 

SMZ 

= 

Scale factor for Z-coordi nates 

The 

following card 

is 

repeated NLOADS times N=1,NL0ADS 

Card 12 (2F12.6) 



Column 



Contents 

1-12 

AL0ADS{1,N) 

= 

Load increment for first load case 

13-24 

AL0ADS(2,N) 

= 

Initial load state before applying increment 

Card 13 is omitted 

if 

lELET = 0 

Card 13 (F12.6) 



Column 



Contents 

1-12 

DELTOT 

= 

Constant temperature for non-thermal loading 

Card 14 is omitted 

if 

lELEM = 0 

Card 14 (F12.6) 



Column 



Contents 

1-12 

DELMOT 

= 

Constant moisture content for non-hygro 




scopic loading 

The following cards 

are repeated NDIFM times 

K=1 

,NDIFM (Cards 15-28) 

Card 15 

(5E12.6) 



Column 



Contents 

12 

EK11(K,1) 

= 

E-j-j tension modulus 

13-24 

EK11(K,2) 

= 

E-J 1 compression modulus 

25-36 

EK22{K,1) 

IZ 

^22 modulus 

37-48 

EK22(K,2) 

= 

Egg compression modulus 

49-60 

EK33(K,1) 

= 

Egg tension modulus 

61-72 

EK33(K,2) 

s 

Egg compression modulus 
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Card 16 

3E12.6 



Column 


Contents 


1-12 

GK23(K) 

= G 23 modulus 


13-24 

GK13(K) 

= G-J 2 modulus 


25-36 

GK12(K) 

= G ^2 modulus 


Card 17 

(6E12.6) 



Column 


Contents 


1-12 

SP1(K,1) 

= Elastic limit stress for - 

E-j tension 

13-24 

N1(1,K,1) 

= Ramberg-Osgood coefficient n-j 

for a-j - e-j 



tension 


25-36 

K 1 ( 1 ,K, 1 ) 

= Ramberg-Osgood coefficient 

for a-j - E-j 



tension 


37-48 

SPI1(K,1) 

- Bilinear intersect stress for 

a-j - e-j 



tension 

49-60 

N1(2,K,1) 

= Ramberg-Osgood coefficient 

for a-j - £-, 



tension 


61-72 

K 1 ( 2 ,K, 1 ) 

= Ramberg-Osgood coefficient K 2 

for a-j - e-j 



tension 


Card 18 (6E12.6) 



Column 


Contents 


1-12 

SP1(K,2) 

= 


13-24 

N1(1,K,2) 

= 


25-36 

K 1 ( 1 ,K, 2 ) 

= Same as Card 17 but for - 

^1 

37-48 

SP1(K,2) 

= compression 

49-60 

N1(2,K.2) 

= 


61-72 

K 1 ( 2 ,K, 2 ) 
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Card 19 (6E12.6) 



Column 



Contents 

1-12 

SP2{K,1) 

- 


13-24 

N2(1.K,1) 

= 


25-36 

K2(1,K,1) 

=5 

Same as Card 17 but for Og 

37-48 

SPI2(K,1) 

S 

tension 

49-60 

N2{2,K,1) 

= 


61-72 

K2(2,K.l) 

= 


Card 20 (6E12.6) 



Column 



Contents 

1-12 

SP2(K,2) 

= 


13-24 

N2(1,K,2) 

8 


25-36 

K2(1,K,2) 

S 

Same as Card 17 but for 

37-48 

SPI2(K,2) 

S 

compression 

49-60 

N2(2,K,2) 

= 


61-72 

K2(2,K,2) 



Card 21 

(6E12.6) 



Column 



Contents 

1-12 

SP33(K,1) 

= 


13-24 

N33(1,K,1) 

= 


25-36 

K33(1,K,1) 

- 

Same as Card 17 but for ■ 

37-48 

SPI33(K,1) 

= 

tension 

49-60 

N33{2,k,l) 

= 


61-72 

K33(2,K,1) 

=: 







Contents 


Card 22 (6E12.6) 
Column 


1-12 

SP33(K,2) 

= 

13-24 

N33(1,K,2) 


25-36 

K33(1,K,2) 

= Same as Card 17 but for 

37-48 

SPI33(K,2) 

= compression 

49-60 

N33(2,K,2) 


61-72 K33(2,K,2) 

Card 23 (6E12.61 
Column 

Contents 


1-12 SP23(K) 

13-24 N23(1,K) 

25-36 K23(1,K) = Same as Card 17 but for T 23 - Y 23 

37-48 SPI23(K) = shear 

49-60 N23(2,K) 

61-72 K23(2,K) 

Card 24 (6E12.6) 


Column 

1-12 

SP13(K) 

Contents 

13-24 

N13(1,K) 

= 

25-36 

K13(1,K) 

= Same as Card 17 but for 

37-48 

SPI13(K) 

= shear 

49-60 

N13(2,K) 

= 

61-72 

K13(2,K) 
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Card 25 (6E12.6) 




Column 



Contents 


1-12 

SP3(K) 

s 



13-24 

M3(1,K) 

- 



25-36 

K3(1,K) 

= Same as Card 17 but for " Y-ig 

37-48 

SPI3(K) 

= shear 



49-60 

N3(2,K) 




61-72 

K3(2,K) 

= 



Card 26 (5E12.6) 




Column 



Contents 


1-12 

SLIO.K) 

= Ultimate 

stress for a-j 

- tension 

13-24 

SL1(2,K) 

= Ultimate 

stress for cj-j 

- e-j compression 

25-36 

S120,K) 

= Ultimate 

stress for ^2 

- Eg tension 

37-48 

SL2(2,K) 

= Ultimate 

stress for Og 

- Eg compression 

49-60 

SL33(1,K) 

= Ultimate 

stress for 

- £2 tension 

61-72 

SL33(2,K) 

= Ultimate 

stress for 

- E^ compression 

Card 27 

(3E12.6) 




Column 



Contents 


1-12 

SL23(K) 

= Ultimate 

stress for T 23 

■ ^23 

13-24 

SL13(K) 

= Ultimate 

stress for x^^ 

- ^13 

25-36 

SL3(1,K) 

= Ultimate 

stress for x-j 2 

- ^12 

Card 28 (6E12.6) 




Column 



Contents 


1-12 

UK12(K.l) 

= Poisson's 

ratio v ^2 

tension 

13-24 

UK12{K,2) 

= Poisson's 

ratio v ^2 

compression 


25-36 

UK23(K,1) 

= Poisson's ratio in tension 

37-48 

UK23(K,2) 

= Poisson's ratio ^23 in compression 

49-60 

UK13(K,1) 

= Poisson's ratio v-j^ in tension 

61-72 

UK13(K,2) 

= Poisson's ratio v-j^ in compression 

If 

no thermal analysis is required skip to card 40 

The 

following cards 

are repeated NDIFM times 

K =1 

,NDIFM (Cards 29-39) 

Card 29 (6112) 


Column 


Contents 

1-12 

NTl(K) 

= Number of linear segmented points for 
modulus percent retention curve 

13-24 

NT2(K) 

= Number of linear segmented points for 
Egg modulus percent retention curve 

25-36 

NT33(K) 

= Number of linear segmented points for 
Egg modulus oercent retention curve 

37-48 

NT23(K) 

= Number of linear segmented points for 
Ggg modulus percent retention curve 

49-60 

NT13(K) 

= Number of linear segmented points for 
G-jg modulus percent retention curve 

61-72 

NT3(K) 

= Number of linear segmented points for 
G-J 2 modulus percent retention curve 

Card 30 (3112) 


Column 


Contents 

1-12 

NT4(K) 

= Number of linear segmented points for 


a-j thermal coefficient curve 
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= Number of linear segmented points for 
ct 2 thermal coefficient curve 
= Number of linear segmented points for 
thermal coefficient curve 

Card 31 (6F12.0) I=UNT1(K) 

Column Contents 

1-12 PERMRl(l 3 K) = Percent retention of modulus 

at point I 

13-24 TEMPI ( I, K) = Temperature at point I 

etc. Repeated NTl(K) times 
Card 32 (6F12.0) I=1,NT2(K) 

Col unin Contents 

1-12 PERMR2(I,K) = Same as card 31 but for E 22 
13-24 TEMP2{r,K) = modulus 

Card 33 (6F12.0) 1=1 ,NT33(v) 


13-24 NT5(K) 

25-36 NTALP3(K) 


Column Contents 

1-12 PMR33(I,K) = Same as card 31 but for E^^ 

13-24 TMR33(I,K) = modulus 

Card 34 (6F12.0) I=1,NT23(K) 


^ol unin Contents 

1-12 PMR23(I,K) = Same as card 31 but for 

13-24 TRM23(I,K) = modulus 

Card 35 (6F12.0) I=1,NT13(K) 


Column 

1-12 


Contents 


13-24 


PMRi2(I,K) 

TMR13(I,K) 


Same as card 31 but for Gno 

13 

modulus 
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^^M(6F12.0) W.nT3(K) 
Column 


1-12 

PERMR3(I,K) 

= Same as card 31 but for G, 


TEMP3(I,K) 

12 

~ modulus 

Card 37 3(E12.S,F12.o) 

I=1,NT4(K) 

Column 





Contents 

1-12 

ALP1(I,K) 

thermal coefficient at point 

13-24 

TEMP4(I,K) 

= Temperature at point I 

etc. 

Repeated NT4(K) 

times 

Card 38 

3(E12.5, F12.0) 

I=1,NT5(K) 

Col umn 





Contents 

1-12 

ALP2{I,K) 

Same as card 37 but for 

13-24 

TEHP5{I,K) 

coefficient 

Card 39 3(E12.5.R2.0) 

I=1,NTALP3(K) 

Column 





Contents 

1-12 

flLP3(I,K) 

Same as card 37 but for a 

13-24 

TflLP3(I,K) 

3 

coefficient 


" "" s«p t. „ 

The feHowing cards ar® repeated NDIFM tJires 
K^I^NDIFM (Cards 40-50) 

Card 40 (6112) 


Column 

1-12 

13-24 


NM1(K) 


NM2(K) 


Contents 

- Number of linear segmented points ■ 
En modulus percent retention curve 

- Number of linear segmented points f 


T 6 . 




modulus percent retention curve 


25-36 

NM33(K) 

= Number of linear segmented points for 
modulus percent retention curve 

37-48 

NM23(K) 

= Number of linear segmented points for 
G 22 modulus percent retention curve 

49-60 

NM13(K) 

Number of linear segmented points for 
G-J 3 modulus percent retention curve 

61-72 

NM3(K) 

= Number of linear segmented points for 
^12 percent retention curve 

Card 41 

. (3112) 


Column 


Contents 

1-12 

NM4(K) 

= Number of linear segmented points for 
e-j coefficient 

13-24 

NM5(K) 

= Number of linear segmented points for 
Bg coefficient 

25-26 

NBETA3(K) 

= Number of linear segmented points for 
B 3 coefficient 

Card 42 

(6F12.0) 1=1 

,NM1(K) 

Column 


Contents 

1-12 

PERMR4(I,K) 

= Percent retention of modulus 

at point I 

13-24 

TEMM1(I,K) 

= Moisture content at point I 

etc. 

Repeated NMl(K) times 
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Card 43 (6F12.0) I=1,NM2{K) 


Column Contents 

1-12 PERMR5(IjK) = Same as card 42 but for E22 

13-24 TEMM2(I,K) = modulus 

Card 44 (6F12.0) I=1,NM33(K) 

Col unin Contents 

1-12 PMMR33(I,K) = Same as card 42 but for 

13-24 TMMR33(I,K) = modulus 

Card 45 (6F12.0) I=1,NM23(K) 

Column ‘ Contents 

1-12 PMMR23(I,K) = Same as card 42 but for G23 

13-24 TMMR23(I,K) = modulus 

Card 46 (6F12.0) I=1,NM13(K) 

Column Contents 

1-12 PMMR13{I,K) = Same as card 42 but for G-j^ 

13-24 TMMR13(IsK) = modulus 

Card 47 (6F12.0) I=1,NH3(K) 

1-12 PERMR6(I,K) = Same as card 42 but for G-jg 

13-24 TEMM3(I,K) = modulus 

Card 48 3(E12.S, F12,0) I=1,NH4(K) 


Col umn Contents 

1-12 BETA1(I,K) = hygroscopic coefficient at 

point I 


13-24 TEMM4(I,K) = Moisture content at point I 


Repeated NM4(K) times 


etc. 
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Card 49 3(E12.5. Fl?.n1 

I=1,NM5(K) 

Column 

Contents 

1-12 BETA2(I,K) 

= Same as card 48 but for 

13-24 TEMM5(I,K) 

- coefficient 

Card 50 3fE12.5.F12.n) 

I=1,NBETA3(K) 

Column 

Contents 

1-12 BETA3(I,K) 

= Same as card 48 but for 

13-24 TBETA3(I,K) 

= coefficient 

Card 51 f6F12.6l 

Column 

Contents 

1-12 THE(l) 

= Angle number 1 in degrees 

13-24 THE(2) 

= Angle number 2 in degrees 

etc. Repeated NANG times 

The following card is repeated NDS times 

1=1, NDS 

Card 52 (4I3.2F12.nl 

Column 

Contents 

1-3 INODED(I) 

= I 

4-6 IN0DE{I,1) 

= U - displacement code 

7-9 IN0DE(I,2) 

= V - displacement code 

10-12 IN0DE(I,3) 

= W - displacement code 
= 1 for force or non-zero displacement 


boundary condition 

= 2 for prescribed zero-displacement 


I 
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13-24 YY(I) 

Y coordinate of node I before being scaled 


by SMY 

25-36 ZZ(I) 

Z coordinate of node I before being scaled 


by SMZ 

The following card is repeated NE times 

1=1 ,NE 


Card 53 (6X,5I6l 


Column 

Contents 

1-6 

= Blank 

7-12 ND(I,1) 

= Node number 1 of element I 

13-18 ND(I,2) 

= Node number 2 of element I 

19-24 ND(I,3) 

= Node number 3 of element I 

25-30 IMAT(I) 

= Material number of element I 

31-36 ITHETA(I) 

= Angle number of element I 

Card 54 121121 


Column 

Contents 

1-12 NDCST 

Number of non-zero displacement constraints 

13-24 NFCST 

Number of non-zero force constraints 

If NDCST = 0 skip 

to card 56 

The following card 

is repeated NDCST times 

1=1, NDCST 


Card 55 (2I12.F12.n) 


Column 

Contents 

1-12 NODED(I) 

= Node number of constrained node 

13-24 MODE 

= Code for constraint 
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-r 




= 2 for V constraint 
= 3 for W constraint 

25-36 DCST{I) = Displacement constraint increment 

If NFCST = 0 no more input is required 
The following card is repeated NFCST times 
1=1, NFCST 

Card 56 (2I12,F12.0) 


Column 

1-12 NODEF(I) 

13-24 MODE 

25-36 FCST(I) 


Contents 

Same as card 55 but for force 
constraints 


Notes: 


GENERAL: Input units need only be consistent except for thermal 

properties which must be degrees F. 

Card 6 NE<400, NDS<400, NDIFM<10, NANG<10 

Card 7 NL0ADS< 5 

Card 10 The first and last load increments are always printed 

Card 11 The scale factors are multiplied by the Y and Z coordinates 

to obtain the final Y and Z coordinates 
Card 12 AL0ADS(2,N) is applicable to thermal and hygroscopic 

loading only, AL0ADS(1,N) must be input as 0.0 for inplane 
loadings. 

Card 32 Node numbers must be given in counter-clockwise order 


A listing of the computer program is available upon request. 
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